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H.  W.  Alt 

Universitdt  Bonn,  Germany 

Entropy  in  non-equilibrium  thermodynamics 

We  give  an  answer  to  the  question  whether  the  entropy  principle  in  rational  thermodynamics 
can  be  derived  from  more  basic  principles.  We  formulate  an  ’’evolution  principle”  and  show  that  in 
certain  classical  situations  it  is  equivalent  to  the  entropy  principle.  We  then  discuss  cases  related 
to  order  parameters  and  interfaces. 


Alexis  Bonnet 

Paris,  France 

On  the  number  op  singular  points  on  a  free-boundary 

On  two  particular  problems,  I  will  describe  how  one  can  get  informations  on  the  number  of 
pieces  and  singular  points  of  a  free-boundary.  The  first  problem  comes  from  the  minimization  of  the 
Mumford-Shah  functional  in  computer  vision.  The  second  set  of  results,  which  have  been  derived 
by  Regis  Monneau  (CERMICS,  PARIS),  deals  with  the  obstacle  problem  in  two  dimensions. 


Andrejs  Cebers 

Institute  of  Physics,  Latvia 

2D  Free  surface  flows  of  magnetic  liquids 

Competing  long-range  magnetic  repulsion  and  surface  tension  forces  are  causing  the  devel¬ 
opment  of  remarkable  labyrinthine  patterns  in  Hele-Shaw  cells.  The  control  parameter  (mag¬ 
netic  Bond  number)  of  those  pattern  formation  phenomena  is  determined  by  the  ratio  of  the 
self-magnetic  field  and  capillary  forces.  In  the  frame  of  the  simple  model  accounting  for  the  self- 
magnetic  field  energy  of  magnetic  liquid  in  Hele-Shaw  cell  and  surface  energy  of  free  boundaries 
the  threshold  values  of  the  magnetic  Bond  number  for  various  instability  modes  in  different  ge¬ 
ometries  (circular  droplet,  infinite  stripe  etc.)  have  been  calculated.  The  patterns  developing  at 
supercritical  values  of  the  magnetic  Bond  number  are  determined  by  the  numerical  solution  of 
the  free  boundary  problem  for  the  magnetic  liquid  motion  in  Hele-Shaw  cell  under  the  action  of 
magnetic  and  capillary  forces.  Here  the  corresponding  numerical  results  obtained  for  the  magnetic 
liquid  Hele-Shaw  flows  by  applying  boundary  integral  equation  technique  are  given.  Numerical 
simulation  of  the  pattern  formation  in  simplest  case  of  enough  small  magnetic  Bond  numbers 
illustrate  the  establishment  of  nonsymmetric  figures  of  equilibrium  (dumbbell,  three  lobe  figure). 
It  is  shown  by  numerical  simulation  results  that  development  of  the  figures  of  equilibrium  with 
the  symmetry  axis  of  the  higher  order  than  three  is  impossible  due  to  vertex  splitting  instability. 
No  peculiarities  in  the  energy  time  dependence  curves  are  noticed  at  vertex  splitting  events.  It  is 
shown  by  numerical  simulation  that  after  several  vertex  splitting  events  at  intermediate  stages  the 
long- living  transient  configuration  close  to  Steiner  trees  are  formed.  Under  longer  time  periods 
evolution  of  those  configurations  to  the  simpler  tree-like  structures  takes  place.  That  evolution 
opposite  to  the  first  stage  is  driven  by  the  surface  tension  forces.  The  influence  of  the  viscos¬ 
ity  contreist  on  those  phenomena  is  studied  by  the  comparison  with  the  case  of  the  labyrinthine 
instability  development  of  the  bubble  in  infinite  plane  layer  of  magnetic  liquid. 

As  it  is  possible  to  conclude  from  the  numerical  simulation  results  of  the  labyrinthine  instability 
of  the  magnetic  liquid  drops  for  the  patterns  formed  stripes  of  the  constant  width  are  characteristic. 
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From  the  analysis  of  the  properties  of  the  stripe  structures  several  very  interesting  conclusions  can 
be  drawn: 

1)  in  the  case  of  the  infinite  system  of  the  straight  magnetic  stripes  in  the  Hele-Shaw  cell  effective 
surface  tension  of  the  equilibrium  pattern  due  to  the  long-range  magnetic  forces  is  exactly  equal 
to  zero; 

2)  at  enough  fast  magnetic  field  increase  undulation  instability  of  the  stripe  system  leading  to  the 
formation  of  chevron  patterns  exists; 

3)  the  similar  concept  of  the  negative  effective  surface  tension  can  be  applied  as  well  to  the 
description  of  the  disordered  systems  of  the  magnetic  liquid  stripes,  for  example,  magnetic  liquid 
foams; 

4)  various  defects  in  the  magnetic  liquid  stripe  systems  can  be  described  on  the  basis  of  2D  smectic 
approach. 

Several  theoretical  conclusions  mentioned  above  are  illustrated  by  the  numerical  simulation  of  the 
corresponding  free  boundary  problems.  So  for  magnetic  liquid  foam  structure  at  the  magnetic 
Bond  number  increase  the  undulation  instability  of  the  single  stripes  is  observed.  It  is  illustrated 
also  by  the  numerical  simulation  of  the  single  stripe  that  the  development  of  the  overextensions 
leading  to  the  rupture  of  the  stripe  is  impossible  due  to  the  vertex  splitting  instability  as  the  result 
the  formation  of  the  alternating  finger  patterns. is  observed. 

Finally  at  the  end  of  the  first  part  several  further  issues  of  the  labyrinthine  pattern  formation 
are  considered; 

a)  topological  instabilities  at  the  dumbbell  relaxation  at  switching  off  the  magnetic  field  and 
displacement  of  the  stripe; 

b)  peculiarities  of  the  influence  of  the  long-range  magnetic  interactions  on  Rayleigh- Taylor  and 
Saffman- Taylor  instabilities. 

In  the  second  part  some  results  for  the  2D  free  boundary  Stokes  flows  of  the  magnetic  liquid 
in  ac  fields  are  given.  As  it  is  known  from  theoretical  investigations  for  2D  magnetic  fluid  droplet 
the  instability  with  respect  to  the  formation  of  the  star  fishes  exists.  By  numerical  simulation 
of  the  corresponding  free  boundary  problem  it  is  shown  that  2-spike  shape  has  the  lowest  energy 
and  is  usually  obtained  at  the  relaxation  from  arbitrary  initial  shape.  3-  spike  shape  can  exist 
as  metastable  state.  Shapes  with  higher  number  of  the  spikes  are  impossible  due  to  the  vertex 
splitting  instability.  According  to  the  numerical  simulation  results  different  regimes  of  the  magnetic 
liquid  droplet  motion  are  possible  in  low  frequency  rotating  magnetic  fields.  At  low  frequencies  the 
droplet  is  elongated  and  follows  synchronously  the  field.  Behaviour  at  higher  frequencies  depends 
on  the  value  of  magnetic  Bond  number  -  for  magnetic  Bond  number  greater  as  critical  one  at 
high  frequencies  complex  regimes  of  the  droplet  motion  including  forward  and  backward  rotations 
arise.  Simple  analytical  model  for  the  droplet  behaviour  in  rotating  field  can  be  derived  which 
shows  good  agreement  with  the  numerical  simulation  data.  With  that  model  devils  staircase  for 
the  droplet  rotation  frequency  in  dependence  on  elliptically  polarized  field  rotation  frequency  is 
obtained. 

Joint  work  with  I.Drikis  and  S.Lacis. 


Richard  Ewing 

Texas  A  &  M  University,  USA 

Free  boundary  problems  arising  in  applications  of  fluid  flow 

IN  POROUS  MEDIA 

Understanding  the  fate  and  transport  of  contaminants  to  determine  water  quality  and  to  de¬ 
velop  remediation  strategies  or  optimizing  the  recovery  of  hydrocarbons  in  petroleum  applications 
each  require  the  ability  to  model  multiphase  flow  in  heterogeneous  three-dimensional  reservoirs. 
Model  equations  and  corresponding  parameters  must  be  determined  at  the  appropriate  length 
scales  to  describe  the  scaled  physics  of  flow.  The  interaction  of  the  various  phases  in  multiphase 
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flows  yields  very  sharp  fluid  interfaces  that  travel  accross  the  reservoir  under  the  influence  of  in¬ 
jection  and  production  rates.  The  speed  of  these  fronts  depends  in  a  highly  nonlinear  way  on  the 
dynamic  state  of  the  fluids  and  fluid  pressures.  These  free  boundary  problems  must  be  treated 
computationally  due  to  the  strong  nonlinearities  and  couplings.  Effective  simulators  require  accu¬ 
rate  numerical  methods  on  general  geometries.  Use  of  mixed  finite  element  methods  and  local  grid 
refinement  will  be  discussed.  Example  calculations  for  field  simulations  in  aquifers  or  reservoirs 
with  complex  boundaries  will  be  presented.  Parallelization  of  the  codes  will  also  be  discussed. 


Antonio  Fasano 

Universita  di  Firenze,  Italy 

Two  SCALE  DIFFUSION  PROBLEMS 

In  many  problems  of  practical  interest  diffusion  of  mass  and/or  heat  occurs  in  a  medium  having 
a  fine  component  in  which  the  same  process  takes  place,  possibly  accompained  by  other  phenom¬ 
enon,  usually  coupled  with  the  bulk  process.  A  remarkable  example  is  a  heat  conducting  medium 
undergoing  nucleation  of  crystals,  in  which  the  whole  crystallization  process  is  crucially  influenced 
by  temperature  and  is  accompained  by  the  release  of  latent  heat.  Another  interesting  case  is  the 
freezing  of  droplets  of  liquid  dispensed  in  another  liquid  with  lower  melting  temperature.  The 
droplets  phase  change  effects  the  heat  diffusion  process  on  the  large  scale  and  in  turn  the  solidifi¬ 
cation  of  each  droplet  is  governed  by  the  local  value  of  the  bulk  temperature.  Filtration  through 
porous  media  can  also  provide  peculiar  examples  of  two-scale  diffusion  with  free  boundaries,  if 
the  solid  matrix  contain  ’’active”  components  (e.g.  swelling  moisture  absorbing  grains).  An  ex¬ 
tremely  complex  process  exhibiting  two-scale  diffusion  is  polymerization  of  gaseous  monomers  in 
an  agglomerate  of  ting  catalytic  particles.  Here  we  have  to  consider  the  monomer  diffusion  in  the 
agglomerate  (the  large  scale),  and  through  the  polymer  shell  growing  around  the  particles  (the 
small  scale),  the  latter  being  a  free  boundary  problem.  Not  only  it  is  important  to  consider  the 
simultaneous  heat  diffusion  (polymerization  is  highy  exothermic),  but  also  the  expansion  velocity 
fields  in  the  two  scales  must  be  reconstructed.  A  review  of  mathematical  models  will  be  presented. 


Carlos  Kenig 

University  of  Chicago,  USA 

MONOTONICITY  formulas  and  REGULARITY  OF  TWO  PHASE 
FREE  BOUNDARY  PROBLEMS 

We  will  discuss  the  role  of  monotonicity  formulas  in  the  proof  of  regularity  of  free  boundaries 
in  two  phase  situations.  We  will  emphasize  some  recent  works,  on  parabolic  problems  (joint  with 
L.Caffarelli),and  on  elliptic  problems  of  Poisson  type  (joint  with  L.Caffarelli  and  D.Jerison).The 
results  have  applications  to  combustion  theory  (in  the  first  case)  and  to  fluid  flows  (in  the  second 
case). 
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Masayasu  Mimura 

The  Univeristy  of  Tokyo,  Japan 


During  the  past  decade,  various  interesting  phenomena  of  patterns  and  interfaces  have  been 
observed  in  chemical  reaction  processes,  solidifications,  combustion,  binary  alloy,  liquid  crystals  in 
the  fields  of  natural  sciences  and  engineering.  In  order  to  understand  the  dynamics  of  patterns, 
most  of  them  are  described  by  certain  types  of  free  boundary  problems.  Along  this  line,  I  would  like 
to  demonstrate  several  biological  pattern  formation  arising  segregation  or  aggregation  of  biological 
individuals,  which  are  described  by  interfacial  or  free  boundary  problems.  The  tool  used  there  is 
a  singular  limit  analysis. 


Stanley  Osher 

University  of  California  Los  Angeles,  USA 

The  level  set  method  for  analyzing  and  computing 

FREE  BOUNDARY  PROBLEMS 

In  1987,  together  with  J.  A.  Sethian,  we  devised  a  new  numerical  procedure  for  capturing 
the  motion  of  fronts.  The  method  uses  a  fixed  grid  and  finds  the  front  as  a  particular  level  set, 
moving  with  time,  of  a  scalar  function.  The  technique  has  since  been  applied  to  an  enormous 
number  of  topics,  ranging  from  computer  vision  to  computational  fluid  dynamics  and  beyond. 
The  method  handles  topological  merging  and  breaking,  works  in  any  number  of  space  dimensions, 
does  not  require  that  the  moving  front  be  written  as  a  function,  captures  steep  gradients  and  cusps 
in  the  front  (using  modern  capturing  schemes,  e.g.,  ENO)  and  is  easy  to  program.  Theoretical 
justification,  involving  the  notion  of  viscosity  solutions,  has  been  given.  Many  applications  and 
extensions  have  been  recently  developed.  For  this  talk  we  will  describe  the  method  applied  to 
multiphase  compressible  interface  motion,  Stefan  problems,  and  problems  in  materials  science. 
Additionally,  the  level  set  point  of  view  enables  us  to  obtain  some  new  analytic  results  concerning 
crystal  growth  and  shape. 


Constantine  Pozrikidis 

University  of  California  at  San  Diego,  USA 

Dynamic  simulations  of  the  flow  of  suspensions  of 

LIQUID  DROPS 

Suspensions  of  liquid  drops  are  encountered  in  a  variety  of  chemical  engineering  and  bioengi¬ 
neering  systems,  and  have  been  studied  using  a  host  of  theoretical  and  experimental  methods. 
In  this  research,  we  carry  out  dynamic  simulations  of  the  flow  of  suspensions  of  two-dimensional 
liquid  drops  with  constant  interfacial  tension  in  simple  shear  flow.  The  computations  rely  on  the 
method  of  interfacial  dynamics,  which  is  an  advanced  implementation  of  the  Boundary  Integral 
Method.  The  numerical  task  is  the  iterative  solution  of  a  system  of  second  kind  integral  equations. 
The  convergence  is  expedited  by  deflating  the  spectrum  of  the  double-layer  operator.  The  results 
allow  us  to  characterize  the  motion  in  terms  of  effective  rheological  properties,  and  also  study 
the  dynamics  of  the  evolving  microstructure.  The  pair  distribution  function  for  suspensions  of 
moderate-viscosity  drops  is  shown  to  be  significantly  different  from  that  of  suspensions  of  rigid 
particles,  and  this  is  explained  in  terms  of  the  ability  of  the  drops  to  slide  over  each  other  during 
collision.  The  seemingly  random  motion  of  one  individual  drop  in  the  suspension  is  described 
in  terms  of  an  effective  hydrodynamic  self-diffusivity  which  is  computed  from  the  results  of  the 
simulations. 
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Mike  Savage 

University  of  Leeds,  UK 

Meniscus  roll  coating  and  the  bead  break  instability 

Free  boundary  problems  are  ubiquitous  in  liquid  film  coating.  In  this  presentation  attention 
is  focused  on  Meniscus  Roll  Coating  -  a  coating  regime  which  arises  when  the  inlet  is  starved  so 
that  flow  rate  is  small  and  a  ’’fluid  bead”  is  located  in  the  nip  bounded  by  two  moving  rolls  and 
two  free  surfaces.  The  pressure  field  is  entirely  subambient  dominated  by  capillary  pressure  at  the 
upstream  meniscus.  The  flow  consists  of  a  fluid  transfer  jet  or  ’’snake”  circumnavigating  two  large, 
closed  recirculations  in  order  to  transport  inlet  flux  to  the  upper  roll.  The  stability  of  the  two  - 
dimensional  base  flow  to  small  amplitude  perturbations  is  considered.  It  is  shown  that  ,  as  upper 
roll  speed  is  gradually  increased,  the  upstream  meniscus  passes  through  the  nip  and  instability  in 
the  form  of  bead  break  soon  follows. 


Jose  Alexandre  Scheinkman 

University  of  Chicago,  USA 

Some  free  bounday  problems  that  arise  in  economics  and  finance 

In  this  lecture  we  will  discuss  two  examples  of  free  boundary  problems  arising  in  economics.  The 
first  example  arises  in  the  theory  of  asset  pricing  in  financial  markets.  The  second  example  deals 
with  dynamic  optimization  in  the  presence  of  irreversible  decisions,  and  is  inspired  by  questions 
on  the  economics  of  environmental  management. 


H.A.  Stone 

Harvard  University,  USA 

Free-boundary  problems  in  fluid  dynamics  with  emphasis 

ON  VISCOUS  FLOWS 

An  overview  of  the  wide  variety  of  free-boundary  problems  that  arise  in  fluid  dynamics  will 
be  given.  Flows  for  which  integral  equation  methods  (for  either  high  or  low  Reynolds  numbers), 
numerical  approaches  valid  for  arbitrary  Reynolds  numbers,  one-dimensional  approximations,  and 
similarity  solutions  are  useful  will  be  summarized  and  their  physical  applications  indicated.  Also, 
cases  where  usually  distinct  analytical  and  numerical  approaches  are  combined  (e.g.  integral 
equation  methods  with  the  lubrication  approximation)  will  be  demonstrated.  Several  examples 
drawn  from  our  own  research  (and  also,  of  course,  the  research  of  others)  in  viscous,  multiphase 
flows  will  illustrate  the  interplay  of  numerics  and  experiments  for  obtaining  additional  physical 
insight,  including  coalescence  of  deformable  drops  in  dilute  suspensions,  drop  and  thread  breakup, 
and  the  conical  shape  of  the  interface  between  two  dielectric  liquids  exposed  to  an  electric  field. 
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Claude-Michel  Brauner 

Universite  Bordeaux  I,  France 

A  SURVEY  OF  RECENT  MATHEMATICAL  RESULTS  ON  FBPS  IN  COMBUSTION 
I  will  report  on  the  following  papers; 

1)  L.A  Caffarelli  (Courant  Institute),  C.  Lederman  (Buenos  Aires)  and  N.  Wolanski  (Buenos 
Aires):  Pointwise  and  viscosity  solutions  for  the  limit  of  a  two-phase  parabolic  singular  per¬ 
turbation  problem 

2)  C.  Lederman  and  N.  Wolanski  (Buenos  Aires):  Viscosity  solutions  and  regularity  of  the  free 
boundary  for  the  limit  of  an  elliptic  two-phase  singular  perturbation  problem 

3)  C.  Lederman  (Buenos  Aires),  J.  L  Vazquez  (Univ.  Autonoma  de  Madrid)  and  N.  Wolanski 
(Buenos  Aires):  Uniqueness  of  solution  to  a  free  boundary  problem  from  combustion 

These  papers  deal  with  models  of  propagation  of  premixed  flames  of  the  following  types, 

Au^  —  uf  =  Pe{u^)  or  Au^  =  j3e{u^), 

where  1/e  stands  for  the  large  activation  energy  and  is  a  normalized  temperature,  (de{s)  — 
Convergence  results  have  been  obtained  as  e  -t  0  to  FBPs 

Aw  —  wt  =  0  or  Aw  =  0, 

with  free  boundary  conditions  such  that  u  —  0,  Uu  —  %/2M  (one-phase  problem)  or  (wj)^  — (w[7)^  = 
2M  (two-phase  problem)  where  M  is  related  to  /3. 

In  the  two-phase  parabolic  case,  limits  are  shown  to  be  solutions  to  the  FBP  in  pointwise 
and  viscosity  senses.  In  the  elliptic  case,  among  other  results,  two-phase  nondegenerate  limits  are 
shown  to  be  classical  solutions.  In  the  one-phase  parabolic  case,  different  notions  of  solution  to 
the  FBP  are  proved  to  coincide  and  produce  a  unique  solution. 

4)  A.  Bonnet  (ENS  Paris)  and  F.  Hamel  (Universite  Paris  VI):  Existence  of  non  planar  solutions 
of  a  simple  model  of  premixed  Bunsen  flames 

In  a  Bunsen  burner  with  long  shaped  aperture  the  flame  front  has  a  bidimensional  conical 
shape.  Far  away  from  the  symmetry  axis  it  behaves  like  a  planar  front  with  constant  speed  cq.  It 
is  proved  the  existence  of  a  solution  to 

Au  —  cDyU f{u)  =  0  inR^, 

with  two  different  limits  as  the  radii  go  to  infinity:  w  0  as  r  oo  in  an  open  cone  of  aperture 
angle  a  determined  by  c  =  co/sina,  and  w  — >■  1  in  the  open  complementary  cone. 

5)  J.  Audounet  (Universite  Toulouse  3),  V.  Giovangigli  (Ecole  Polytechnique)  and  J.-M.  Roque- 
joffre  (Universite  Toulouse  3):  An  integral  equation  modelling  the  propagation  of  a  point  source 
initiated  flame 

The  time  evolution  of  the  radius  of  a  lean  spherical  flame,  is  described  by  the  integro- 
differential  equation  (Joulin,  1985) 

Rdi/2B  =  RlogR -i- Eq{t),  i?(0)  =  0, 

where  di/2  denotes  the  time  derivative  of  order  1/2.  The  function  q{t)  represents  the  point  source, 
and  the  parameter  E  represents  its  strength.  It  has  been  proved  that  the  flame  radius  either  tends 
to  0,  or  to  1,  or  to  -foo,  according  to  whether  the  parameter  E  is  below,  equal  to,  or  above  a 
critical  parameter  Ecr{q)-  Asymptotics  of  the  different  regimes  have  also  been  obtained. 

6)  G.S.  Namah  (University  Bordeaux  1)  and  J.-M.  Roquejofffe  (Universite  Toulouse  3):  A  class 
of  nonlinear  parabolic  and  Hamilton- Jacobi  equations  arising  in  solid  propellant  combustion 
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A  qualitative  model  for  the  propagation  of  a  flame  front  in  a  solid  propellant  reads  as  follows. 


ut-Au  =  J?(X,u)^l  +  |Dul2,  X£R^. 

It  is  proved  that  periodicity  is  space  generates  periodic  in  time  solutions.  Related  Hamilton-Jacobi 
equations  are  also  studied  in  the  class  of  viscosity  solutions.  The  corresponding  homogenization 
problems  (by  letting  the  spatial  period  go  to  zero)  have  also  been  analyzed,  to  which  effective 
ode’s  are  derived.  For  this  end,  profit  is  taken  of  the  knowledge  of  the  time-asymptotic  behaviour 
of  the  above  solutions. 

7)  C.M.  Brauner  (Universite  Bordeaux  1),  A.  Lunardi  (Universita  di  Parma)  and  Cl.  Schmidt- 
Laine  (ENS  Lyon):  Stability  and  bifurcation  of  non  planar  travelling  fronts  in  a  two-phase 
problem 

8)  C.M.  Brauner  (Universite  Bordeaux  1),  J.  Hulshof  (University  Leiden)  and  Cl.  Schmidt-Laine 
(ENS  Lyon):  The  saddle  point  property  for  focusing  selfsimilar  solutions  in  a  one  phase  prob¬ 
lem 

In  the  limit  of  high  activation  energy,  certain  models  consist  of  one-  or  two-phase  problems 
such  as 

ut  =  Au-\-  f{u,Du) 

on  both  sides  of  a  free  boundary  a;  =  ^{y,t),  where  u  =  u*  is  given  and  jump  conditions  on  the 
normal  derivative  of  u  are  prescribed.  An  adequate  splitting  of  the  unknown  yields  the  elimination 
of  the  front  and  leads  to  a  fully  nonlinear  evolution  equation.  The  latter  formulation  is  specially 
efficient  in  the  nonlinear  stability  analysis  of  travelling  wave  (TW)  solutions. 

In  a  model  problem  in  dimension  2,  the  planar  TW  becomes  unstable  at  a  critical  value  ul  of 
u*,  at  which  an  infinite  number  of  bifurcated  branches  of  nonplanar  TWs  accumulates. 

Similarly,  the  free  surface  can  be  eliminated  in  the  context  of  spherical  flames  for  studying 
stable  and  unstable  manifolds  of  focusing  self-similar  solutions. 


Amable  Linan 

Universidad  Politecnica  de  Madrid,  Spain 
FBP  IN  DIFFUSION  CONTROLLED  COMBUSTION 
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Bill  Bold 

UMIST,  U.K. 

Non-monotonic  curvature-dependent  propagation  of  a  front 

If  an  interface  propagates  at  a  speed  S{k)  that  depends  non-monotonically  on  the  mean  cur¬ 
vature  K  then  the  shape  of  the  interface  experiences  an  anti-diffusive  instability  whenever  the 
derivative  S'{k)  is  negative — taking  k  to  be  positive  when  concave  in  the  direction  of  propaga¬ 
tion.  While  negative  diffusion  has  connotations  of  non-physical  behaviour,  such  as  violating  the 
second-law  of  thermodynamics,  ill-posedness,  etc.,  it  turns  out  that  if  S'{k)  is  negative  only  in 
bounded  windows  of  curvature-space  then  any  growth  that  arises  from  the  instability  is  nonlinearly 
bounded.  Having  a  “diffusion  coefficient”  that  changes  sign,  as  k  changes,  makes  the  problem  more 
tractable  than  a  problem  with  uniformly  negative  diffusion. 

A  great  deal  can  be  deduced  about  the  appearance  and  behaviour  of  the  instability,  includ¬ 
ing  the  surprising  result  that,  under  very  simple  conditions,  an  interface  can  propagate  with  a 
negatively  diffusive  instability  for  an  infinite  time  and  remain  well-posed  with  respect  to  initial 
conditions!  The  effect  of  noise  is  another  matter  altogether,  and  any  noise  that  decays  more  slowly 
than  a  Gaussian  spectrum  will  cause  an  instant  anti-diffusive  growth  as  soon  as  S'{k)  becomes 
negative  (that  is,  as  normal  propagation  of  the  interface  causes  the  value  of  k  to  move  through 
some  threshold  value).  The  full  growth  of  the  instability  is  still  limited  within  the  windows  of 
anti-diffusive  curvature. 

Wide  classes  of  steadily  propagating  solutions  can  also  be  found  with  piecewise  continuous 
curvature.  These  resemble  a  phase-separation  in  curvature  space  into  positively  and  negatively 
curved  solutions.  For  appropriate  models  of  S{k),  normal  propagation  has  the  effect  of  driving 
solutions  into  the  anti-diffusive  regime  where,  seeded  by  noise,  such  phase  separation  would  take 
place. 

A  simple  generalisation  of  the  Kuramoto-Sivashinsky  equation  should  lead  to  very  similar 
dynamics,  with  the  added  presence  of  a  regularising  term  involving  a  fourth-derivative,  or  Laplacian 
of  curvature.  It  is  not  clear  at  this  stage  how  important  the  details  of  any  such  regularisation  are 
likely  to  be  in  determining  the  dynamics.  Unlike  models  for  separation  into  physical  phases  (as 
opposed  to  phases  of  curvature)  it  is  likely  that  regularisation  would  play  a  secondary  role  which 
does  not  change  the  overall  dynamics  very  much. 

An  interesting  and  rich  dynamic  can  be  expected  for  fronts  which  propagate  with  fairly  simple 
models  for  S{k).  Models  of  this  nature  have  been  predicted  using  slowly,  but  not  too  slowly, 
varying  flame  analyses  of  mixtures  near  stoichiometry,  reacting  via  one-step  Arrhenius  kinetics  at 
large  activation  temperature  with  certain  combinations  of  non-unit  Lewis  numbers.  Calculations 
using  fuller  models  of  flames  at  not  too  weak  curvatures  or  stretch-rates  are  needed  to  identify  if 
the  phenomenon  can  be  found  more  generally.  Because  of  the  instability,  experiments  in  regimes 
of  negative  S'{k)  are  likely  to  become  problematic. 

It  is  worth  pointing  out  that  not  only  flames  can  propagate  at  non-monotonic  speeds  S{k). 
Other  forms  of  interface  possess  this  kind  of  law,  including  the  meandering  of  rivers.  Studying  the 
dynamics  of  this  new  class  of  problem  promises  to  reveal  features  that  will  help  in  understanding 
the  unsteady  behaviour  of  a  range  of  interfaces  in  which  anti-diffusive  instabilities  are  restricted 
to  windows  of  solution  space. 
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R.  Luciano 

Rome  Univ.  II,  Italy 


Pavel  Krejci 

Weierstrass  Institute  for  Applied  Analysis 
and  Stochastics,  Germany 
Hysteresis  in  thermoplasticity 

We  propose  a  model  for  the  behavior  of  elastoplastic  materials,  where  the  material  character¬ 
istics  (elasticity  modulus,  yield  surface)  depend  on  temperature.  The  construction  is  based  on  the 
idea  of  Prandtl  (1928)  and  Ishlinskii  (1944)  to  describe  real  elastoplastic  materials  by  a  parallel 
combination  of  simple  elastic  -  perfectly  plastic  elements  parametrized  by  the  yield  thresholds 
with  a  distribution  function  which  is  to  be  identified  from  experiments.  Rather  than  large  systems 
of  evolution  variational  inequalities,  it  is  more  convenient  to  work  with  their  solution  operators 
which  belong  to  the  class  of  hysteresis  operators.  In  our  thermoelastoplasticity  model,  we  simply 
assume  that  the  Prandtl-lshlinskii  distribution  function  depends  also  on  temperature.  Since  the 
constitutive  law  is  given  in  terms  of  hysteresis  operators,  it  is  natural  to  expect  that  also  the  inter¬ 
nal  energy  and  entropy  ensuring  the  thermodynamical  consistency  of  the  model  will  have  the  form 
of  hysteresis  operators.  Our  main  goal  is  to  prove  the  well-posedness  of  the  corresponding  system 
of  the  equation  of  motion  and  energy  balance  equation.  While  the  modeling  part  is  independent  of 
the  space  dimension,  existence  and  uniqueness  results  for  systems  of  partial  differential  equations 
involving  temperature-dependent  hysteresis  operators  are  obtained  only  in  one  space  dimension 
due  to  restrictions  imposed  by  the  Sobolev  embedding  theorems. 

Joint  work  with  J.  Sprekels. 


Martin  Brokate 

University  of  Kiel,  Germany 

Nonlinear  kinematic  hardening 

We  discuss  constitutive  laws  in  elastoplasticity  of  kinematic  nonlinear  hardening  type;  ’nonlin¬ 
ear’  here  means  that  the  direction  of  the  movement  of  the  yield  surface  differs  from  the  direction 
of  the  plastic  strain  increment  by  a  nonlinear  correction  term.  A  representative  model  goes  back 
to  Armstrong  and  Frederick;  its  multisurface  version,  the  Chaboche  model,  is  widely  used  in  en¬ 
gineering.  We  present  results  concerning  the  wellposedness  and  asymptotic  behaviour  of  such 
models. 
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Vasilios  Alexiades 

University  of  Tennessee,  USA 

Alloy  solidification  with  convection  in  the  melt 

Most  models  of  alloy  solidification  are  severely  limited  by  the  assumption  of  constant  density, 
thus  excluding  all  convective  effects.  We  develop  a  thermodynamicaly  consistent  model  for  binary 
alloy  solidification  that  incorporates  energy,  species  and  momentum  conservation,  constitutional 
supercooling,  as  well  as  temperature  and  concentration  dependence  of  thermophysical  parameters. 
The  crucial  aspect  is  the  development  of  an  equation  of  state  capturing  the  thermochemistry  of 
the  phases.  A  numerical  scheme  will  also  be  outlined. 


Harald  Haller 

Technische  Universitat  Munchen,  Germany 

Micromechanical  models  and  homogenization  of  an  alastic  plate 

CONTAINING  FIBERS  OF  SHAPE  MEMORY  ALLOYS 

Based  on  a  micromechanical  approach  and  on  a  generalized  Achenbach-Miiller  model  for  the 
shape  memory  alloys  a  thermomechanical  model  of  an  elastic  thin  plate  containing  fibers  of  shape 
memory  alloys  is  derived.  Existence  and  uniqueness  of  the  solutions  is  proven. 

Since  the  physical  parameters  (such  as  heat  conductivity,  elasticity  coefficients,  ...)  are  discon¬ 
tinuous  and  oscillate  between  different  values  characterizing  each  of  the  components,  the  numerical 
solution  of  the  derived  equations  requires  a  lot  of  computation  time  and  for  many  applications  it 
is  even  not  possible  to  solve  the  equations  numerically. 

Using  homogenization  techniques  one  can  get  a  good  approximation  of  the  macroscopic  be¬ 
haviour  of  such  heterogeneous  materials  by  letting  the  parameter  e,  which  describes  the  fineness 
of  the  microscopic  structure,  tend  to  zero  in  the  given  equations  and  therefore  the  computation 
time  can  be  reduced.  We  derive  the  homogenized  (or  effective)  equations  for  the  fiber-reinforced 
plate  and  prove  the  convergence  of  the  solutions. 

Numerical  simulations  of  the  amplitude  control  of  the  vibrating  plate  show  a  very  good  corre¬ 
spondence  to  experimental  results. 


Tomas  Roubicek 

Charles  University,  Czech  Republic 

Twinning  and  inelastic  response  of  crystalic 

MATERIALS 

Plasticity  in  crystalic  materials  can  be  caused  beside  slip  also  by  twinning.  This  is  basically 
caused  by  a  dissipative  martensitic  phase-change  created  by  thermal  or  dislocation  mechanisms. 
These  mechanisms  are  activated  if  the  stress  exceeds  certain  limit  so  that  the  potential  barrier 
between  particular  phases  became  partly  ’’transparent”.  In  the  scalar  case,  this  problem  seems  to 
be  modelled  by  the  author’s  microstructure  evolution  model  modified  by  considering  the  dissipation 
energy  nonconstant,  being  activated  by  a  sufficiently  large  stress.  The  phenomena  like  hysteresis, 
observed  experimentally  in  shape-memory  alloys  exposed  to  time-varying  loading,  can  thus  be 
modelled  at  least  for  a  simply  twinned  configuration. 
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Barbara  Wagner 

Technische  Universitdt  Munchen,  Germany 
Numerical  and  asymptotic  results  on  the  linear  stability  of  a 

THIN  FILM  SPREADING  DOWN  A  SLOPE  OF  SMALL  INCLINATION 

For  a  thin  liquid  film  moving  down  a  slope,  the  lubrication  approximation  with  a  small  Navier- 
Slip  is  used  as  a  model.  The  travelling  wave  solution  is  derived  for  small  inclination  angle  a, 
using  singular  perturbation  methods,  and  compared  to  the  numerical  solution.  For  the  linear 
stability  analysis  we  combine  numerical  methods  with  the  long-wave  approximation  and  find  a 
small  but  finite  critical  a*  below  which  the  flow  remains  linearly  stable.  This  is  contrasted  with 
the  vanishing  of  the  hump  of  the  travelling  wave  solution.  Finally,  this  prevailing  of  linear  stability 
of  the  travelling  wave  at  small  inclination  angles  is  compared  with  recent  related  results  using  a 
precursor  model  instead.  Here  though,  a  strong  dependency  on  the  magnitude  of  the  contact  angle 
is  found,  which,  we  think,  has  not  been  observed  before. 
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Daniel  D.  Joseph 

University  of  Minnesota,  USA 

Cavitation  and  the  state  of  stress  in  aflowing  liquid 

The  prob.lem  of  the  inception  of  cavitation  is  formulated  in  terms  of  a  comparison  of  the  breaking 
strength  or  cavitation  threshold  at  each  point  of  liquid  sample  with  principal  stresses  there.  A 
cavity  will  open  in  the  direction  of  the  maximum  tensile  stresses  which  is  45°  from  the  plane  of 
shearing  in  pure  shear  of  a  Newtonian  fluid.  A  criterion  of  maximum  tension  is  proposed  which 
unifies  the  theory  of  cavitation,  the  theory  of  maximum  tensile  strengh  of  liquid  filaments  and  the 
theory  of  fracture  of  amorphous  solids.  Experiments  which  support  these  ideas  are  discussed  and 
some  new  experiments  are  proposed. 


John  Lowengrub 

University  of  Minnesota,  USA 

Cahn-Hilliard  hydrodynamics  and  topological  transitions 

The  competition  between  surface  tension  and  instability  often  leads  to  singularities  in  interfacial 
flows  in  which  interfaces  collide  and/or  self- intersect  and  the  topology  of  the  flow  changes.  Before 
the  singularity,  sharp  interface  models  provide  good  models  of  the  flow.  At  the  singularity  time, 
however,  these  models  break  down  and  additional  features  of  the  flow  must  be  considered  such 
as  the  partial  miscibility  between  “immiscible”  fluid  components.  In  this  talk,  I  will  introduce  a 
new  model  (Cahn-Hilliard  Hydrodynamics)  in  which  sharp  interfaces  are  replaced  by  microscopic 
mixing  layers  that  account  for  the  partial  miscibility  that  real  fluids  always  display.  Analysis  and 
simulations  will  demonstrate  the  model’s  ability  to  capture  changes  in  topology. 


Constantine  Pozrikidis 

University  of  California  at  San  Diego,  USA 

Numerical  studies  of  singularity  formation  at  free  surfaces  and 

FLUID  INTERFACES  IN  TWO-DIMENSIONAL  STOKES  FLOW 

We  consider  the  analytic  structure  of  interfaces  in  several  families  of  steady  and  unsteady  two- 
dimensional  Stokes  flows,  focusing  on  the  formation  of  corners  and  cusps.  Previous  experimental 
and  theoretical  studies  have  suggested  that,  without  surface  tension,  the  interfaces  spontaneously 
develop  such  singular  points.  We  investigate  whether  and  how  corners  ad  cusps  actually  develop  in 
a  time-dependent  flow,  and  assess  the  stability  of  stationary  cusped  shapes  predicted  by  previous 
authors.  The  motion  of  the  interfaces  is  computed  with  high  resolution  using  a  boundary  integral 
method  for  three  families  of  flows.  In  the  case  of  a  bubble  that  is  subjected  to  the  family  of 
straining  flows  devised  by  Antanovskii,  we  find  that  a  stationary  cusped  shape  is  not  likely  to 
occur  as  the  asymptotic  limit  of  a  transient  deformation.  Instead,  the  pointed  ends  of  the  bubble 
disintegrate  in  a  process  that  is  reminiscent  of  tip  streaming.  In  the  case  of  the  flow  due  to  an 
array  of  point-source  dipoles  immersed  beneath  a  free  surface,  which  is  the  periodic  version  of  a 
flow  proposed  y  Jeong  &  Moffatt,  we  find  evidence  that  a  cusped  shape  indeed  arises  as  the  result 
of  a  transient  deformation.  In  the  third  part  of  the  numerical  study,  we  show  that,  under  certain 
conditions,  the  free  surface  of  a  liquid  film  that  is  leveling  under  the  action  of  gravity  on  horizontal 
or  slightly  inclined  surface  develops  an  evolving  corner  or  cusp.  In  certain  cases,  the  film  engulfs 
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a  small  air  bubble  of  ambient  fluid  to  obtain  a  composite  shape.  The  structure  of  a  corner  or  a 
cusp  in  an  unsteady  flow  does  not  have  an  unique  shape,  as  it  does  at  steady  state.  In  all  cases,  a 
small  amount  of  surface  tension  is  able  to  prevent  the  formation  of  a  singularity,  but  replacing  the 
inviscid  gas  with  a  viscous  liquid  does  not  have  a  smoothing  effect.  The  ability  of  the  thin-film 
lubrication  equation  to  produce  mathematical  singularities  at  the  free  surface  of  a  leveling  film  is 
also  discussed. 


Leonid  Antanovskii 

Monash  University,  Australia 

Hysteresial  behaviour  of  a  pointed  drop  in  Taylor’s 

FOUR-ROLLER  MILL 

Time-dependent  deformation  of  a  two-dimensional  inviscid  drop  subjected  to  quasi-steady 
straining  flow,  such  as  that  created  in  Taylor’s  four-roll  mill  [Taylor,  1934]  at  slowly  varying 
rotation  rate  of  the  cylinders,  is  addressed.  Following  the  spirit  of  matched  asymptotic  expan¬ 
sions,  the  problem  is  decomposed  into  the  outer  Stokes  problem  with  no  drop  effect,  which  is  solved 
by  the  boundary-element  method,  and  the  inner  free-boundary  problem  for  the  drop  behaviour  in 
unbounded  fluid.  Using  complex  variable  technique,  a  broad  class  of  explicit  solutions  is  found, 
which  are  described  by  a  rational  conformal  mapping  of  the  unit  disc  onto  the  flow  domain,  with 
time-dependent  coefflcients  satisfying  ordinary  differential  equations.  Some  numerical  simulations 
of  the  resulting  equations  are  implemented,  which  suggest  that,  with  increasing  rotation  rate  of 
the  cylinders,  a  sudden  cusping  of  a  sufficiently  small  drop  occurs  at  a  critical  strain.  This  cusped 
configuration  survives  up  to  any  strain,  but  when  the  strain  is  being  gradually  removed,  the  drop 
becomes  suddenly  rounded  at  a  lower  critical  strain.  This  proves  that,  under  certain  circumstances, 
the  behaviour  of  a  2D  inviscid  drop  placed  within  a  quasi-steady  creeping  flow  is  hysteresial. 
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Ralf  Kornhuber 

Universitdt  Stuttgart,  Germany 

On  the  fast  solution  of  free  boundary  problems 

A  wide  range  of  free  boundary  problems  occurring  in  engineering  and  industry  can  be  rewritten 
as  a  minimization  problem  for  a  strictly  convex,  piecewise  smooth  but  non-dilferentiable  energy 
functional.  The  algebraic  solution  of  the  related  discretized  problem  is  a  very  delicate  question, 
because  usual  Newton  techniques  can  not  be  applied. 

We  propose  a  new  approach  based  on  convex  minimization.  Each  iteration  consists  of  one  step 
of  a  globally  convergent  descent  method  and  an  additional  energy  reducing  correction.  While  the 
initial  (fine  grid)  smoothing  provides  global  convergence  of  the  overall  iteration,  the  subsequent 
(coarse  grid)  correction  is  intended  to  accelerate  the  convergence  speed.  We  present  a  general  con¬ 
vergence  theory  involving  asymptotic  multigrid  convergence  rates  and  discuss  several  applications. 


Maurizio  Paolini 

Universitd  di  Udine,  Italy 

Approximation  of  crystalline  motion  by  mean  curvature 

BY  A  NONLINEAR  REACTION- DIFFUSION  EQUATION 

We  present  some  results  concerning  the  so-called  crystalline  mean  curvature  flow  approximated 
by  a  reaction-diffusion  equation  which  formally  reads  as  foil  ows: 

e^dtu  -  e^divT(Vu)  -f  f{u)  —  eg{x,t) 

where  T  is  a  maximal  monotone  graph,  homogeneous  of  degree  1,  describing  the  underlying  crys¬ 
talline  anisotropy;  /  is  the  derivative  of  a  bistable  potential,  and  p  is  a  given  forcing  term.  A 
few  numerical  simulations  suggest  good  approximation  properties  similar  to  those  obtained  for  the 
classical  isotropic  mean  curvature  flow. 
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Michael  J.  Shelley 

New  York  University,  USA 

Interface  instability  and  dynamics  in  non-Newtonian  Hele-Shaw  flow 

The  flow  of  shear-thinning  liquids  in  thin  gaps  is  important  to  such  applications  as  injection 
molding  (polymeric  liquids),  and  display  device  design  (liquid  crystals).  From  the  point  of  view 
of  pattern  formation,  the  displacement  of  a  shear-thinning  liquid  by  gas  in  a  Hele-Shaw  cell  can 
give  interfacial  patterns  -  with  the  suppression  of  tip-splitting,  and  appearance  of  side-branching 
-  that  are  very  reminiscent  of  anisotropic  solidification  systems,  though  no  explicit  anisotropy 
is  present.  To  study  such  problems,  we  have  developed  the  natural  generalization  of  Darcy’s 
law  for  non-Newtonian  liquids,  beginning  with  fluid  mechanical  models  such  as  the  generalized 
Navier-Stokes  equations,  and  an  Oldroyd  fluid  model,  both  of  which  allow  for  shear-thinning, 
and  a  non-monotonic  stress/strain  relation.  Unlike  the  Newtonian  Darcy’s  law,  the  viscosity  is 
now  dependent  upon  the  squared  pressure  gradient,  and  leads  to  a  nonlinear  elliptic  boundary 
value  problem  for  the  fluid  pressure.  We  are  studying,  both  analytically  and  numerically,  the 
modification  of  the  Saffman- Taylor  instability  of  an  expanding  gas  bubble  in  a  shear-thinning 
liquid.  We  have  shown  that  shear-thinning  leads  to  enhanced  length  scale  selection.  We  are  also 
developed  numerical  methods  for  simulating  the  dynamics  of  the  gas/liquid  interface,  and  our 
simulations  show  the  nonlinear  development  of  patterns  composed  of  stabilized  fingers  that  do  not 
tip-split. 

This  is  joint  work  with  Ljubinko  Kondic  (Courant)  and  Peter  Palffy-Muhoray  (Liquid 
Crystal  Institute). 


Alfred  Schmidt 

Univ.  Freiburg,  Germany 

A  POSTERIORI  ERROR  ESTIMATION  AND  ADAPTIVE  METHODS  FOR 

STEFAN  PROBLEM 

We  present  an  adaptive  finite  element  method  with  a  posteriori  error  estimation  for  degenerate 
parabolic  problems  like  the  two-phase  Stefan  problem 

dtu  -t-  A/3(u)  =  f. 

The  a  posteriori  error  estimates  are  based  on  duality  techniques,  involving  degenerate  para¬ 
bolic  problems  in  nondivergence  form.  They  give  computable  bounds  for  the  enthalpy  error  in 
L°°(0,T;  ff“^(f2))  and  for  the  temperature  error  in  L‘^{0,T-,  L^{Cl)).  The  adaptive  method  tries  to 
generate  quasi-optimal  meshes  by  equidistribution  of  the  local  error  estimates  in  space  and  time. 
Numerical  results  in  two  and  three  space  dimensions  will  be  presented. 

(joint  work  with  R.H.  Nochetto  and  C.  Verdi) 
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Pierre  Gremaud 

North  Carolina  State  University,  USA 

New  results  in  numerical  conservation  laws 

In  this  talk,  new  results  related  to  the  numerical  approximation  of  scalar  conservation  laws 
will  be  presented.  In  particular,  the  role  played  not  only  by  the  size  of  the  grid,  but  also  its 
structure,  i.e.,  uniform  vs.  nonuniform,  Cartesian  vs.  nonCartesian,  will  be  carefully  analysed. 
The  underlying  question  of  the  proper  notion  of  consistency  will  also  be  discussed.  Various  new 
results  for  specific  numerical  methods  will  be  presented. 

The  approach  has,  among  other  things,  led  to  an  explanation  of  the  supraconvergence  phe¬ 
nomenon,  and  to  optimal  error  estimates  in  the  truly  multidimensional  case.  To  the  authors’ 
knowledge,  both  types  of  results  are  unique  in  the  framework  of  conservation  laws. 

Joint  work  with  B.  Cockburn  and  X.  Yang. 
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P.H.  Gaskell 

The  University  of  Leeds,  UK 

Stability  analysis  for  reverse  roll  coating 

Two  mathematical  models  are  presented  for  meniscus  reverse  roll  coating  in  which  there  is 
steady  flow  of  a  Newtonian  fluid  in  the  narrow  gap,  or  nip,  between  two  co-rotating  rolls  in  the 
absence  of  body  forces. 

A  model  based  on  lubrication  theory  and  incorporating  visco-capillary  boundary  conditions 
gives  predictions  for  the  pressure  field,  the  meniscus  locations  and  flow  rate  as  a  function  of 
roll  speed  ratio,  capillary  number  and  a  geometry  parameter.  Solutions  are  found  to  be  highly 
dependent  on  roll  speed  ratio  and  degree  of  inlet-starvation. 

A  second  model  incorporates  the  nonlinear  free-surface  boundary  conditions  and  the  presence 
of  a  dynamic  wetting  line,  adjacent  to  the  web  on  the  upper  roll.  The  latter  requires  the  imposition 
of  an  apparent  contact  angle  and  a  slip  length  so  that  numerical  solutions  for  the  velocity  and 
pressure  fields  over  the  entire  domain  can  be  obtained  using  a  finite  element  formulation  of  the 
governing  equations  of  motion.  Predictions  via  both  models,  for  film  thicknesses  attached  to  the 
web  and  lower  roll  as  functions  of  roll  speed  ratio,  are  in  close  agreement  with  each  other  and  also 
with  corresponding  experimental  data. 

Finite  element  solutions  reveal  that  the  flow  domain  (fluid  bead)  consists  of  large  recirculations 
and  fluid  transfer-jets,  and  that  the  flow  structure  is  continuously  transformed  as  the  roll  speed 
ratio  and  inlet  flux  are  varied.  For  a  given  capillary  number  and  geometry  (roll  radius  to  gap  height) 
two  critical  flow  rates  can  be  identified  in  the  starved/  ultra-starved  regimes,  each  dependent  on 
roll  speed  ratio,  and  corresponding  to  the  appearance/  disappearance  of:  (i)  a  secondary  transfer- 
jet;  (ii)  a  saddle  point  at  the  nip  with  a  downstream  sub-eddy.  Predicted  flow  structures  and 
critical  flow  rates  are  found  to  be  in  qualitative  agreement  with  experimental  observation. 

Finally,  it  is  shown  that,  for  a  given  inlet  flux,  there  is  a  maximum  roll  speed  ratio  such  that 
only  below  this  value  can  two-dimensional  steady-state  solutions  be  obtained  via  the  finite  element 
method.  It  is  suggested  that  this  is  related  to  the  bead-break  instability  which  arises  after  the 
wetting  line  has  passed  through  the  nip.  The  relevant  information  is  display  on  an  associated 
control-space  diagram. 


J.  Eggers 

Universitdt  Gesamthochschule  Essen,  Germany 
Physics  of  drop  formation 

One-dimensional  models  have  proved  to  be  very  useful  to  study  the  breakup  of  liquid  jets. 
They  can  be  derived  from  the  Navier-Stokes  equation  either  by  expanding  the  velocity  field  in 
the  radial  direction  or  by  averaging  over  thin  slices  of  the  fluid.  The  singularity  of  the  equations 
of  motion  which  describes  the  radius  of  the  fluid  jet  going  to  zero  locally  is  characterized  by  a 
balance  of  surface  tension,  viscous,  and  inertial  forces.  The  singular  motion  is  independent  of  fluid 
parameters  and  of  initial  conditions.  Away  from  asymptotics,  one  dimensional  models  are  able  to 
describe  an  even  greater  variety  of  scaling  behavior.  Finally,  the  mechanisms  of  main  drop  and 
satellite  formation  are  studied  in  a  simulation  of  a  stationary,  decaying  jet. 
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K.D.  Danov 

University  of  Sofia,  Bulgaria 

Surfactant  dynamics  in  thin  films 

Surfactants  play  an  important  role  in  many  technological  processes.  They  adsorb  at  the  in¬ 
terfaces  and  reduce  the  interfacial  mobility  and  the  surface  elasticity  due  to  the  change  of  the 
interfacial  tension.  A  hydrodynamic  model  of  the  multi-component  mixture  of  surfactant  includ¬ 
ing  the  description  of  the  bulk-  and  inter-phases  is  discussed.  New  constitutive  equations  for  the 
surface  diffusion  and  heat  fluxes,  shear  and  dilatational  viscosities  are  obtained  using  the  Onsager 
theory  for  the  surface-excess  dissipation  function.  For  that  purpose  the  non-equilibrium  ther¬ 
modynamics  of  a  mixture  of  surfactants  in  the  presence  of  chemical  reactions  between  them  is 
developed.  These  relations  close  the  system  of  equations  describing  the  physical  behavior  of  the 
moving  interface. 

The  mathematical  difficulties  arising  from  the  special  type  of  the  boundary  conditions  are 
illustrated.  The  possibility  of  formation  of  tangential  heat,  diffusion,  and  viscous  boundary  layers 
on  interfaces  close  to  the  three  phase  contact  lines  are  demonstrated.  Non-solved  mathematical 
and  numerical  problems  with  special  technical  and  biological  applications  are  presented. 

The  physical  approach  for  integral  balance  of  mass  and  momentum  in  surface  boundary  layers 
is  used  to  simplify  the  system  of  equations.  The  applicability  of  this  approach  is  illustrated  on 
the  computation  of  the  meniscus  profile  in  coating  system.  The  influence  of  surface  elasticity, 
viscosity,  dynamic  and  static  contact  angles,  slide  velocity,  and  pressure  vacuum  on  the  position 
of  the  free  three  phase  contact  line  is  investigated. 
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Bjorn  Gustafsson 

Royal  Institute  of  Technology,  Sweden 

An  exponential  trnsform  and  regularity  of  free  boundaries 

IN  TWO  dimensions 


I  will  discuss  a  new  tool  for  proving  regularity  of  the  free  boundary  for  obstacle  type  problems 
with  analytic  data  in  two  dimensions.  Although  this  is  perhaps  the  simplest  nontrivial  case  of  a 
regularity  problem  for  free  boundaries,  a  complete  solution  was  not  given  until  1991,  by  M.  Sakai. 
For  the  corresponding  problem  in  higher  dimension  the  regularity  problem  is  yet  not  completely 
solved,  as  far  as  I  know.  The  best  known  results  are  due  to  L.  Caffarelli. 

The  new  tool  in  question  is  an  exponential  transform,  which  is  a  function  of  two  complex 
variables  originally  created  by  R.  W.  Carey  and  J.  D.  Pincus  in  an  operator  theoretic  context. 
Precisely  it  is 


E{z,  w)  =  exp[- 


-f  — 

TT  in  (C  - 


dA{0 


TT  Vn  (C  -  z)(C  -  w) 


Outside  the  open  set  0,  E{z,  w)  is  anlytic  in  2:  antianalytic  in  w  and  the  main  result  states  that 
if  n  (locally)  is  the  noncoincidence  set  for  the  obstacle  problem,  or  if  more  generally  the  Cauchy 
transform  of  0  has  an  analytic  extension  from  outside  across  50,  then  E{z,  w)  has  an  analytic- 
antianalytic  extension,  call  it  F{z,w),  across  50.  From  this  one  gets  a  canonical  description  of 
the  free  boundary  50  as  being  the  zero  set  of  the  real  analytic  function  F{z,  z). 

One  naturally  wonders  about  possibility  of  generalizing  the  above  ideas  to  higher  dimensions.  It 
does  not  seem  immediate  how  to  proceed,  but  it  certainly  is  an  interesting  task  for  future  research. 

The  connection  between  the  exponential  transform  and  free  boundary  problems  was  discovered 
rather  recently  by  Mihai  Putinar,  and  all  results  I  will  present  are  obtained  in  collaboration  with 
him. 


Stefan  Luckhaus 

Universitdt  Bonn,  Germany 

Stefan  problems  with  the  gibbs  Thomson  law  are  gradient 
flows  for  the  entropy 

Stefan  problems  in  this  context  are  heat  and  mass  diffusion  systems  with  two  or  more  phases. the 
free  boundary  is  the  phase  interface,  the  gibbs  thomson  law  is  the  additional  interface  condi¬ 
tion, apart  from  the  continuity  equations.it  says  that  either  the  mean  curvature  of  the  interface  or 
the  difference  of  mean  curvature  and  velocity  is  proportional  to  the  jump  in  free  energy  across  the 
interface  divided  by  the  absolute  temperature. this  free  boundary  system  can  formally  be  regarded 
as  a  gradient  flow  for  the  total  entropy.this  leads  to  natural  approximation  schemes, and  allows  in 
addition  to  give  a  stochastic  analogon,a  jump  process  on  a  finite  lattice. 
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Augosto  Visintin 

Universitd  degli  Studi  di  Trento,  Italy 
Nucleation  and  Growth 

Surface  tension  effects  occur  in  several  physical  phenomena;  in  particular,  they  are  responsible 
for  the  high  undercooling  required  for  solid  nucleation  in  a  purely  liquid  system.  (For  several 
metals  this  undercooling  may  even  reach  the  order  of  hundreds  of  degrees,  see  e.g.  [1,2].) 

In  a  solid- liquid  system  these  effects  can  be  represented  by  the  classical  Gibbs-Thomson  law, 
which  prescribes  the  relative  temperature  9  to  be  proportional  to  the  mean  curvature  k  of  the 
interface  S  between  phases; 

9  —  —CK  on  5,  (1) 

see  e.g.  [3].  Due  to  the  smallness  of  the  coefficient  c,  this  requires  the  use  of  a  mesoscopic  length 
scale  (typically  of  the  order  of  10“^  cm),  see  e.g.  [5]. 

Evolution  is  then  represented  by  the  law  of  mean  curvature  flow  with  a  forcing  term: 

av  =  9  +  CK  on  5  (2) 


where  v  is  the  normal  velocity  of  the  interface,  -a  is  a  constant  >  0. 

Here  we  want  to  account  not  only  for  smooth  surface  motion  but  also  for  nucleation  and  other 
discontinuities  in  phase  evolution.  To  that  aim  we  replace  (2)  by  an  equation  of  the  form 


a{v)  =  9  +  CK  on  S  (2) 

where  a  is  bounded,  monotone,  continuous  function  R  R  (with  a'(0)  =  0),  see  [4]. 

A  mesoscopic  model  of  phase  transition  in  two-phase  systems  is  then  obtained  by  coupling  the 
latter  equation  with  the  energy  balance: 


in  the  sense  of  distributions, 


(3) 


where  we  set  x  ~  in  solid,  x  ~  1  in  the  liquid;  Cy  is  the  specific  heat  and  k  the  heat 
conductivity. 

The  resulting  problem  is  nonconvex,  because  of  the  constraint  x  =  ±1-  This  is  at  variance 
with  the  weak  (macroscopic)  formulation  of  the  Stefan  problem,  see  e.g.  [6],  in  which  (3)  is  coupled 
with  the  condition 


X  =  — 1  where  9  <  0,  -1  <  X  <  1  where  0  =  0,  x  =  1  where  0  >  0.  (4) 

Values  of  x  intermediate  between  —  1  and  1  correspond  to  a  mixture  of  liquid  and  solid,  a  so-called 
mushy  region;  this  is  excluded  in  the  previous  (mesoscopic)  model. 

Existence  of  a  solution  has  been  proved  for  the  system  consisting  of  equations  (2)  and  (3) 
coupled  with  appropriate  initial  and  boundary  conditions,  see  [7].  In  general  the  uniqueness  of  the 
solution  is  excluded,  because  of  the  nonconvexity. 
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Mete  Soner 

Carnegie  Mellon  University,  USA 
Option  Pricing  with  Constrains 

I  will  discuss  the  problem  of  determining  the  minimum  cost  of  super-replicating  a  claim  when 
there  are  convex  constraints  on  the  portfolio  weights.  The  resulting  minimal  dominating  price  is 
the  solution  of  a  quasi-variational  inequality.  In  some  interesting  examples,  such  as  the  European 
Call  option  with  short  selling  constraint  or  the  up  and  out  Call  option,  I  will  show  that  the 
minimal  price  is  equal  to  the  price  of  a  related  claim  without  constraints.  This  is  joint  work  with 
M.  Broadie  and  J.  Cvitanic  of  Columbia  University. 


Agnes  Tourin 

Universite  Paris-Dauphine,  France 

Numerical  Approximation  to  some  Variational  Inequalities 
ARISING  IN  Mathematical  Finance 

We  consider  here  several  portfolio  selection  models  leading  to  Variational  Inequalities  and  pro¬ 
pose  second-order  numerical  approximations  to  these  Variational  Inequalities  in  order  to  compute 
the  free  boundaries  between  the  transaction  regions.  The  convergence  of  the  schemes  based  on 
recent  ideas  developed  by  Lions  and  Souganidis  in  the  framework  of  first-order  Hamilton- Jacobi 
equations  is  still  an  open  problem. 

This  work  has  been  done  in  collaboration  with  T.  Zariphopoulou. 


Ruediger  Frey 

ETH  Zurich,  Switzerland 

Optimal  Stopping:  The  Free  Boundary  Approach  and 
Applications  to  Derivative  Asset  Analysis 

Optimal  stopping  problems  are  of  great  importance  in  derivative  asset  analysis,  as  they  are 
closely  related  to  the  pricing  of  American  options.  We  therefore  take  the  problem  of  pricing  the 
American  put  option  as  starting  point  for  our  talk.  We  show  that  the  associated  optimal  stopping 
problem  can  be  formulated  as  a  free  boundary  problem  for  the  value  function  of  the  American  put 
option.  This  relationship  between  free  boundary  problems  and  optimal  stopping  problems  is  of  a 
general  nature.  The  free  boundary  formulation  will  help  us  to  determine  the  Riesz  decomposition 
of  the  value  process  of  the  option  and  hence  the  hedge  portfolio  for  the  American  put.  Finally  we 
sketch  a  new  application  of  optimal  stopping  to  the  superreplication  of  derivatives  in  stochastic 
volatility  models:  It  can  be  shown  that  for  a  large  class  of  European  options  the  value  of  a 
superreplicating  hedge  portfolio  is  given  by  the  price  of  the  corresponding  American  option  in  a 
model  with  constant  volatility  and  hence  by  the  value  of  an  optimal  stopping  problem. 
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BRIEF  INTRODUCTION 

In  the  last  decade  phase-field  models  have  received  a  good  deal  of  attention  from  people  working 
in  continuum  mechanics  and  thermodynamics  as  well  as  from  applied  mathematicians  experts  in 
PDE.  In  particular,  it  is  known  that  phase  field  problems  find  noteworthy  applications  to  different 
phase  transition  phenomena  including  solid-liquid  phase  transitions  and  martensitic  transforma¬ 
tions.  Moreover,  the  asymptotic  analysis  of  related  systems  leads  to  some  interesting  and  partly 
unsolved  free  boundary  problems.  The  aim  of  this  session  is  then  providing  the  state  of  the  art 
and  reviewing  recent  results  addressed  to  these  problems. 
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PHASE  FIELD  MODELS  -  PART  1 


This  first  part  is  focused  on  phase  transition  and  phase  separation  models,  collecting  both 
theoretical  and  numerical  results  on  some  PDE’s  systems  which  extend  the  well-known  Caginalp 
and  Penrose-Fife  approaches.  The  following  talks  will  be  delivered. 


Nobuyuki  Kenmochi 

University  Chiba,  Japan 

Phase  separation  model  with  constraint 

In  his  contribution,  N.  Kenmochi  discusses  phase  separation  models  which  are  derived  from  a 
class  of  free  energy  functionals  of  Ginzburg-Landau  type.  These  models  give  rise  to  a  system  of 
two  parabolic  equations  satified  by  the  (absolute)  temperature  and  the  conserved  order  parameter. 
The  first  equation  comes  from  the  energy  balance  and  is  a  nonlinear  second  order  parabolic  PDE 
including  singularity  and  at  the  same  time  degeneracy,  while  the  second  is  a  nonlinear  fourth  order 
parabolic  PDE  including  a  maximal  monotone  graph  in  the  expression  of  the  chemichal  potential 
difference.  It  will  be  shown  that  the  main  part  of  this  phase  separation  model  can  be  expressed 
as  the  subdifferential  of  a  lower  semi-continuous  convex  function  defined  on  an  adequate  product 
space,  so  that  the  related  problem  can  be  handled  within  the  abstract  theory  of  evolution  equations 
in  Hilbert  spaces. 


Maurizio  Grasselli 

Politcnico  di  Milano,  Italy 

Asymptotic  analyses  of  a  phase  field  model  with  memory 

M.  Grasselli  deals  with  a  (standard)  phase  field  model  in  which  the  heat  flux  is  given  by  a 
constitutive  assumption  of  Gurtin-Pipkin  type.  The  resulting  nonlinear  system  which  rules  the 
evolution  of  the  relative  temperature  and  of  the  non-conserved  order  parameter  consists  in  a 
hyperbolic  integrodifferential  equation  coupled  with  a  parabolic  variational  inequality.  Several 
questions  related  to  the  well-posedness  of  the  associated  Cauchy-Neumann  problem  have  been 
treated  jointly  with  P.  Colli  and  G.  Gilardi.  Some  asymptotic  analysis  results  are  presented  in  the 
talk,  in  particular;  (i)  convergence  to  a  hyperbolic  phase  relaxation  problem  when  the  interfacial 
energy  coefficient  tends  to  zero;  (ii)  convergence  to  a  quasi-stationary  phase  field  model  when 
the  time  relaxation  parameter  tends  to  zero;  (iii)  convergence  to  a  standard  parabolic  phase  field 
model  when  the  memory  kernel  tends  to  the  Dirac  mass  in  a  suitable  way.  In  all  the  cases  rigorous 
convergence  theorems  are  proved.  Moreover,  in  (ii)  and  (iii),  some  error  estimates  are  also  shown. 
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Cedric  Dupaix 

Universite  de  Franche-Comte, France 

Upper-semicontinuity  of  the  attractorfor  a  singularly 

PERTURBED  PHASE  FIELD  MODEL 

The  contribution  of  C.  Dupaix  is  concerned  with  boundary  value  problems  for  a  phase  field 
model  with  polynomial  nonlinearity  and  for  two  limiting  equations,  namely  the  viscous  Cahn- 
Hilliard  equation  and  the  Cahn-Hilliard  equation.  All  three  problems  possess  a  maximal  attractor. 
The  upper-semicontinuity  of  attractors  is  shown  and  a  lower-semicontinuity  property  is  presented 
for  the  one-dimensional  viscous  Cahn-Hilliard  equation.  These  results  arise  from  a  collaboration 
with  D.  Hilhorst,  I.  Kostin,  and  Ph.  Laurengot.  Morever,  in  the  talk  the  case  of  a  logarithmic 
nonlinearity  (which  is  more  realistic  in  some  contexts  such  as  solid-solid  phase  transitions)  is  con¬ 
sidered  for  the  phase  field  model  and  the  upper-semicontinuity  of  the  global  attractor  is  extended 
to  this  framework. 


Olaf  Klein 

WIAS  -  Berlin,  Germany 

A  PHASE-FIELD  SYSTEM  WITH  SPACE-DEPENDENT  RELAXATION  COEFFICIENT 

O.  Klein  considers  a  phase-field  system  of  Penrose-Fife  type,  with  a  space-dependent  relax¬ 
ation  coefficient.  A  convergent  time-discrete  scheme  is  discussed  for  the  approximation  of  the 
corresponding  initial-boundary  value  problem.  Existence,  uniqueness,  and  convergence  results  are 
derived,  even  for  a  system  with  degenerate  relaxation  coefficient.  In  conclusion,  some  numerical 
methods  and  tests  are  presented. 
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PHASE  FIELD  MODELS  -  PART  2 


The  second  part  is  mainly  devoted  to  the  asymptotic  behaviour  of  PDE’s  systems  with  respect 
to  suitable  kinetic  parameters  and  to  the  characterization  of  the  limit  problems  and  of  their 
solutions.  Indeed,  the  concerned  titles  are  as  follows. 


Danielle  Hilhorst 

Universite  Paris  -  Sud,  Prance 

A  CHEMOTAXIS-GROWTH  MODEL 

D.  Hilhorst  analyses  an  advection-reaction-diffusion  system  proposed  by  Mimura-Tsujikawa 
and  describing  chemotaxis  with  growth.  Chemotaxis  is  a  phenomenon  observed  in  biology,  in 
which  biological  individuals  migrate  towards  higher  gradients  of  some  chemical  substance  -  the 
chemotactic  substance  -  which  is  produced  by  themselves.  The  two  unknown  functions  are  the 
density  of  biological  individuals  and  the  concentration  of  chemotactic  substance.  As  a  small 
parameter  tends  to  zero,  the  solution  of  the  chemotaxis  problem  converges  to  the  solution  of  a  free 
boundary  problem,  in  which  the  motion  of  the  interface  is  partly  induced  by  its  mean  curvature. 
The  limiting  moving  boundary  problem  possesses  a  unique  smooth  solution  locally  in  time.  Such 
convergence  for  the  solutions  of  the  chemotaxis-growth  problem  has  been  shown  in  a  joint  work 
with  A.  Bonami,  E.  Logak  and  M.  Mimura.  In  addition,  it  has  been  pointed  out  that  properties 
such  as  the  inclusion  and  the  convexity  of  interfaces  are  not  conserved  in  time. 


Elisabeth  Logak 

Ecole  Normale  Superieure,  France 

Mass  conserved  Allen-Cahn  equation  and  area  preserving 

MEAN  CURVATURE  FLOW 

A  mass  conserved  Allen-Cahn  equation  is  considered  by  E.  Logak.  This  equation  has  been 
proposed  by  Rubinstein  and  Sternberg  to  describe  phase  separation  in  binary  mixtures  and  it  can 
be  viewed  as  a  limit  of  the  viscous  Cahn-Hilliard  equation.  The  existence  of  a  unique  smooth 
solution  to  the  area  preserving  mean  curvature  flow  is  assumed  in  some  time  interval.  Under  this 
condition,  we  will  see  that  the  solutions  of  the  mass  conserved  Allen-Cahn  equation  converge  to  this 
limiting  motion  problem,  this  result  following  from  a  cooperation  with  X.  Chen  and  D.  Hilhorst. 


34 


Reiner  Schatzle 

Universitdt  Bonn,  Germany 

Mean  curvature  flow  approximated  by  the  extended 

Schatzle  communicates  on  a  joint  work  with  D.  Hilhorst  and  L.  Peletier  regarding  the  Fisher- 
Kolmogorov  equation  dtu  +  —  Aa  +  e~‘^F'{u)  =  0.  Here,  F{t)  ;=  —  1)^/4  is 

a  double-well  potential  and,  for  7  =  0,  the  above  PDE  reduces  to  the  ordinary  Allen-Cahn 
equation.  Formal  asymptotics  show  that  the  limit  as  e  -4  0  yield  the  mean  curvature  flow.  But 
a  difficulty  arises  since  it  is  not  clear  whether  the  stationary  waves  (that  is,  the  solutions  of  the 
ODE  ')U""  —  U"  +  F'{U)  —  0  with  appropriate  boundary  conditions)  are  unique.  There  could 
also  be  multi-bump  solutions.  Peletier  and  Troy  have  proved  that  monotone  and  odd  solutions  are 
unique  when  0  <  7  <  1/8.  The  main  result  of  this  talk  establishes  that  the  area  functional  is  the 
F— limit  of  £  times  the  stationary  energy  functional  associated  with  the  above  Fisher-Kolmogorov 
equation.  Also,  in  one  space  dimension  the  limit  interfaces  move  with  constant  velocity. 


Giuseppe  Savare 

Istituto  di  Analisi  Numerica  del  CNR,  Italy 

Compactness  properties  for  the  asymptotics  of  some  phase 

field  equations 

The  talk  of  G.  Savare  addresses  the  following  typical  problem,  which  occurs  in  passing  to  the 
limit  in  some  phase  field  models;  for  two  families  of  space-time  dependent  functions  {ue},  {^e} 
(representing,  e.g.,  temperature  and  phase  variable)  we  know  that  the  sum  Ue  +  Ve  converges 
in  some  L^-space  and  that  the  time  integrals  of  a  suitable  space-functionals  evaluated  on  Ue,  Ve 
are  uniformly  bounded  with  respect  to  e.  Can  we  deduce  that  Ue  and  Ve  converge  separately? 
Luckhaus  (1991)  gave  a  positive  answer  to  this  question  in  the  framework  of  the  two-phase  Stefan 
problem  with  Gibbs-Thompson  law  for  the  melting  temperature.  Plotnikov  and  Starovoitov  (1993) 
proposed  a  possible  general  solution  employing  the  original  idea  of  Luckhaus  and  a  space  compact¬ 
ness  assumption.  We  show  how  this  problem  is  related  to  general  metric  and  scaling  properties  of 
the  functionals  and  the  topologies  involved,  providing  sufficient  conditions  for  its  solvability.  This 
approach  allows  simp  le  applications  to  concrete  situations  not  covered  by  the  previous  results, 
like  phase-field  equations  with  Neumann  boundary  conditions  on  the  temperature. 
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DISCUSSION  SESSION 

FREE  BOUNDARY  PROBLEMS  IN  THE  THEORY  OF 

FILTRATION 


ORGANIZED  BY 

S.  Antontsev 

Universidade  da  Beira  Interior,  Portugal 

and 

D.  Hilhorst 

Universite  Paris  -  Sud,  France 


BRIEF  INTRODUCTION 

The  study  of  filtration  problems  is  very  important  both  from  the  theoretical  as  from  the  ap¬ 
plication  points  of  view.  There  are  natural  applications  in  environmental  sciences  when  modelling 
pollution  contaminants  and  their  transport  by  groundwater  and  in  petrol  engineering.  The  math¬ 
ematical  models  which  arise  lead  to  difficult  problems  for  systems  of  nonlinear  partial  differential 
equations  and  free  boundary  problems.  The  scope  of  this  session  is  to  discuss  the  construction 
of  mathematical  models  for  multiphase  filtration.  Such  topics  were  also  discussed  in  the  previous 
congresses  of  this  series  and  are  far  from  being  completely  understood  from  the  mathematical 
point  of  view.  Questions  such  as  existence,  uniqueness  and  qualitative  properties  of  the  solutions 
of  these  free  boundary  models  as  well  as  their  numerical  computation  will  be  raised.  The  speakers 
will  be  the  following. 


LIST  OF  SPEAKERS 

S.N.  Antontsev,  Universidade  da  Beira  Interior,  Portugal 
M.  Chipot,  Universitat  Zurich  -  Irchel,  Switzerland 
M.  Guedda,  Universite  de  Picardie  Jules  Verne,  Prance 
J.  Kacur,  MFF  Comenius  University,  Slovakia 

R.  Kersner,  Hungarian  Academy  of  Sciences,  Hungary 

S.  Shmarev,  Universidad  de  Oviedo,  Spain 
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S.N.  Antontsev 

Universidade  da  Beira  Interior,  Portugal 

Small  parameters  in  the  problems  of  filtration 

OF  TWO  IMMISCIBLE  FLUIDS 

This  talk,  which  describes  joint  work  with  G.  Gagneux,  is  devoted  to  the  investigation  of 
qualitative  properties  of  weak  solutions  of  the  system 

ds 

—  =  (eia(s)Vs)  +  6'(s)VpVs,  (a(s)  >0,  s  6  [0, 1]) 
at 

Qn 

82— =  ik{s)Vp) ,  (ei  >  0, ei ->•  0) 
at 

where  ej  >  0  are  small  parameters.  The  system  models  the  filtration  of  two  immiscible  weakly 
compressible  fluids  [1],  [2],  [3]. 

The  following  qualitative  properties  of  solutions  are  studied  [3]: 
support  properties  of  weak  solutions  (sj  >  0); 
passages  to  the  limit  in  ej  0. 

References 

[1]  Antontsev  S.N.,  Kazhikhov  A.V  and  Monakhov  V.N.  Boundary  Value  Problems  in  Mechanics 
of  Nonhomogeneous  Fluids.  (Russian.).  Nauka.  Sibirsk.  Otdel.  Novosibirsk.  1983.  320  p. 
(English).  Studies  in  Mathematics  and  its  Applications.  22.  North  Holland.  1990.  309  p. 

[2]  Antontsev  S.N.,  Gagneux  G.  Problemes  mathematiques  ouverts  de  la  theorie  de  la  filtration 
diphasique  Laboratoire  de  Mathematiques  Appliquees,  UPRA  1204 

Universite  de  Pau  et  des  Pays  de  I’Adour,  Pau,  France.  Analyse  Non  Lineare  et  Modelisation, 
96/15,  48  p. 

[3]  Gagneux  G,  Monique  Madaune-Tort.  Analyse  mathematique  de  modeles  non  lineaires  de 
I’ingenierie  petroliere.  Mathematiques  &  Applications.  22.  Springer- Verlag  Berlin  Heidelberg, 
1996,  178  p. 


M.  Chipot 

Universitdt  Zurich  -  Irchel,  Switzerland 

On  the  dam  problem  with  leaky  boundary  conditions 

We  would  like  to  present  new  results  -in  particular  of  uniqueness-  for  the  dam  problem  with 
leaky  boundary  conditions.  A  so  called  leaky  boundary  condition  is  a  condition  of  non  linear  flux 
on  the  bottom  of  the  different  reservoirs  considered. 
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M.  Guedda 

Universite  de  Picardie  Jules  Verne,  France 

Blow  up  and  regularity  of  interfaces  for  an  inhomogeneous 

FILTRATION  EQUATION 

The  evolution  of  interfaces  between  fresh  and  salt  groundwater  in  aquifers  in  coastal  areas  can 
be  modelled  by  the  equation 

p{x)ut  =  (u(l  -  u)Ux)x 

where  v,  represents  the  height  of  salt  water.  We  first  give  a  sufficient  condition  on  the  function  p 
for  the  blow  up  of  the  support  of  solutions  of  namely  {{x,t),  0  <  u{x,t)  <  1}.  Similar  results 

have  been  obtained  by  Kamin,  Kersner  and  Galaktionov  about  the  support  of  solutions  of  the 
equation 

p{x)ut  =  {u'^)xx 

We  then  present  regularity  results  for  interfaces  of  solutions  of  This  is  joint  work  with  D.  Hilhorst, 
M.  A.  Peletier  and  S.  I.  Shmarev. 


J.  Kacur 

MFF  Comenius  University,  Slovakia 
Solution  of  degenerate  convection-diffusion  problems 

BY  THE  method  OF  CHARACTERISTICS 

A  new  approximation  scheme  is  proposed  to  solve  a  nonlinear  degenerate  convection-diffusion 
problems.  The  method  is  based  on  relaxation  schemes  developed  by  the  author  in  previous  papers 
and  on  the  method  of  characteristics  initiated  by  O.  Pironneau,  J.  Douglas  and  T.F.  Russel.  The 
mathematical  model  includes  a  nonlinear  degenerate  parabolic  term  and  also  a  nonlinear  convec¬ 
tive  term.  In  the  elliptic  term  the  diffusion  coefficient  depends  on  the  unknown  function.  The 
convergence  of  the  method  is  proved  and  numerical  experiments  supporting  the  theoretical  results 
are  presented.  These  problems  model  the  transport  of  reactive  contaminant  under  Freundlich 
sorption  isotherm. 


R.  Kersner 

Hungarian  Academy  of  Sciences,  Hungary 
Interface  blow-up  in  inhomogeneous  media 

The  model  problem 

p{x)ut  =  {u^)xx  -  cow”, 
u{x,  0)  =  uo{x) 

supported  for 
(Rome). 


will  be  considered,  where  m,n  >  1,  co  nonnegative  and  p{x)  ~  for  large  |a;|. 
The  initial  datum  wo(®)  is  nonnegative  and  has  compact  support. 

We  address  the  following  questions:  when  does  the  solution  remain  compactly 
all  time  and  when  it  becomes  everywhere  positive  after  some  finite  time. 

Two-sided  estimates  for  interfaces  will  also  be  given.  Joint  work  with  A.Tesei 
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S.  Shmarev 

Universidad  de  Oviedo,  Spain 

On  the  regularity  of  interfaces  in  solutions  of  degenerate 
PARABOLIC  equations.  TECHNIQUE  OF  INTERSECTION  COMPARISON 

We  propose  a  new  approach  for  the  study  of  questions  of  regularity  of  solutions  and  interfaces 
for  a  class  of  degenerate  parabolic  equations  with  one  space  variable.  The  method  is  presented  in 
the  framework  of  the  study  of  equations  of  the  form 

(1) 

in  the  range  of  parameters  m  >  1,  0  <  p  <  1,  subject  to  the  condition  m  +  p  —  2  >  0.  It  is  based 
on  the  idea  that  the  local  analysis  of  solutions  of  equation  (1)  near  an  interface  can  be  done  by 
means  of  intersection  comparison  with  the  family  of  travelling  wave  solutions  of  the  equation. 
The  comparison  technique  called  Intersection  Comparison  or  Lap  Number  roughly  says  that  the 
number  of  sign  changes  between  two  solutions  of  equations  like  (1)  is  nonincreasing  in  time.  We 
are  able  to  show  in  this  way  that  the  behaviour  of  the  solutions  near  the  interface,  and  in  particular 
its  local  velocity,  are  determined  in  first  approximation  by  a  travelling  wave  profile. 

•  We  prove  that  for  continuous,  compactly  supported  and  bell-shaped  initial  data  the  solutions 
of  the  Cauchy  problem  for  equation  (1)  have  Lipschitz  continuous  interfaces  in  any  compact  time 
interval  0  <  ti  <  t  <  t2  <  Te-  Here  Te  <  oo  is  the  instant  when  the  solution  extincts. 

•  We  also  obtain  the  formulae  which  express  the  interface  velocity  in  terms  of  the  solution 
profile,  so-called  interface  equations  ,  for  both  expanding  and  receding  waves.  It  is  shown  that  at 
each  instant  t  £  (0,Te)  the  function  is  right  and  left  differentiable. 

The  method  can  be  applied  in  principle  to  quite  general  diffusion-absorption  equations  of 
the  form  with  monotone  nondecreasing  functions  and  other  similar  equations  with  convective  or 
reactive  terms,  or  with  other  types  of  nonlinear  diffusion. 

The  talk  contains  results  of  joint  work  with  V. A. Galaktionov  and  J.L. Vazquez. 
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DISCUSSION  SESSION 

MATHEMATICAL  FREE  BOUNDARY  PROBLEMS 
IN  GLACIOLOGY 


ORGANIZED  BY 


Kolumban  Hutter 

Technische  Hochschule  -  Darmstadt,  Germany 


Luisa  Santos 

Universidade  do  Minho,  Portugal 


and 

Brian  Straughan 

Glasgow  University,  UK 


BRIEF  INTRODUCTION 

We  propose  to  present  a  discussion  in  various  free  boundary  problems  arising  in  the  geophysical 
context  of  glaciology  involving  permafrost  and  its  melt-phase,  ice  sheet  and  ice  shelf  dynamics, 
and  subglacial  till  behaviour.  Emphasis  is  placed  on  demonstrating  the  occurrence  of  free  moving 
boundaries,  internal  and  interface  behaviour,  phase  changes  and,  partly,  mushy  zone  behaviour. 
We  outline  open  mathematical  problems,  review  the  pertinent  physical  literature  and  sketch  pos¬ 
sible  problems  worthy  of  further  mathematical  investigation. 


LIST  OF  SPEAKERS 

J. I.  Diaz,  Univ.  Complutense  de  Madrid,  Spain 

K. H.  Hutter,  Technische  Hochschule  -  Darmstadt,  Germany 
F.  deirisola,  Universita  di  Roma  “La  Sapienza”,  Italy 
J.M.  Urbano,  Universidade  de  Lisboa,  Portugal 
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J.I.  Diaz 

Universidad  Complutense  de  Madrid,  Spain 
On  a  parabolic-hyperbolic  coupled  system  arising  in  glaciology 

The  mechanism  whereby  large  ice  sheets  can  surge  periodically,  (Heinrich  events,  Hudson  bay 
phenomenosj  Canada)  was  recently  studied  by  Fowler  and  Johnson  (1995).  They  proposed  a  two- 
dimensional  ice  sheet  simplified  model  that  includes  basal  ice  sliding  dependent  on  the  basal  water 
pressure,  which  essential  consists  in  the  following  system; 

ht  -  ^{S  +  Q)^h^'^^\hx\^~^hx^  =a  in  Qt  (1) 


Qx-{S  +  Qf  -  7  +  >  0,  Q  >  0  in  Qt 

(Qx  -  (5  -t-  Q)^  -  7  +  Q  =  0  in  Qt 


(2) 


=  {S  +  Q)^ h^\hx\^  in  Qt  (3) 

the  unknown  variables  are  the  ice  depth  h,  the  accumulated  ice  velocity  and  the  basal  water  flow 
Q.  Prescribing  suitable  boundary  conditions  for  h,  Q  and  an  initial  condition  for  h  Fowler  and 
Schiavi  (1997)  solved  numerically  (1),  (2)  and  (3)  using  a  fully  implicit  backwards  finite  difference 
scheme  for  h  and  an  improved  Euler  method  for  Q  and 

The  main  goal  of  this  work  is  to  present  some  results  on  the  existence  of  solutions  for  the  implicit 
backwards  scheme  system  and  some  numerical  experiences  (Fowler  and  Schiavi,  1997)  for  some 
special  data  illustrating  the  surging  mechanism. 

Joint  work  with  E.  Schiavi. 


K.H.  Hutter 

Technische  Hochschule  -  Darmstadt,  Germany 

The  mathematical  foundation  of  ice  sheet  and  ice  shelf  dynamics 

A  physicist’s  view  of  the  pressing  problems:  We  briefly  describe  the  mathematical  formulation 
of  the  boundary  value  problems  that  arise  in  global  ice  sheet  dynamics  and  describe  the  simplest 
models  of  the  ice  sheet/shelf  responses  to  climate  variations.  We  show  in  particular  for  cold  ice 
and  for  polythermal  ice  which  free  boundary  value  problems  the  geophysicist  would  like  to  solve 
and  what  complications  would  have  to  be  overcome,  to  asure  that  the  problems  are  well  posed. 
We  pose  questions  and  ask  mathematicians  for  answers. 


F.  dellTsola 

Universita  di  Roma  “La  Sapienza”,  Italy 

A  FREE  BOUNDARY  VALUE  PROBLEM  FOR  THE  TILL  BELOW  LARGE  ICE  SHEETS: 
WHY  IS  THE  MATHEMATICAL  PROBLEM  PHYSICALLY  RELEVANT? 
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J.M.  Urbano 

Universidade  de  Lisboa,  Portugal 

On  the  mathematical  analysis  of  a  Valley  Glacier  model 

We  study  a  free  boundary  problem  arising  from  theoretical  glaciology: 


ApU^  =  - 
=0 


1  in 
on 


du^ 


=0 


on  r| 


From  the  mathematical  point  of  view  it  consists  of  a  mixed  boundary  value  problem  for  the  p- 
Laplacian  in  a  domain  that  is  not  known  a  priori  and  is  to  be  determined  so  that  a  prescribed 
mass  flux  condition  is  satisfied,  namely  =  Q.  Our  result  establishes  the  existence  of  a  unique 

weak  solution  to  the  problem  by  means  of  the  careful  analysis  of  the  range  of  an  appropriate  real 
function. 

Joint  work  with  J.  F.  Rodrigues. 
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DISCUSSION  SESSION 

DYNAMICS  OF  SURFACE  GROWTH  AND  PHASE  TRANSITION 

KINETICS 


ORGANIZED  BY 


Ivan  Gotz 

Technische  Universitat  Munchen,  Germany 

and 

Boris  Zaltzman 

Ben-Gurion  University  of  Negev,  Israel 


BRIEF  INTRODUCTION 

This  session  is  to  address  the  interface  dynamics  in  phase  change  and  crystal  growth  on  the 
following  two  aspects.  The  first  are  the  kinetic  effects,  either  in  their  proper  thermodynamic 
context  or  as  a  tool  of  regularization  of  ill-posed  problems  {Ch.  Charach,  P.  Fife  and  I.  Gotz). 
The  other  aspect  of  phase  interface  dynamics  to  be  addressed  is  that  of  grain  at  the  surface  of  a 
polycrystalline  material  (A.  Vilenkin  and  B.  Zaltzman).  This  dynamics  is  modelled  by  evolutionary 
fourth  order  equations  in  a  possible  combination  with  additional  ’’kinetic”  equations,  which  have 
so  far  attracted  relatively  little  attention  by  the  FBP  community. 


LIST  OF  SPEAKERS 

A.  Vilenkin,  Hebrew  University  of  Jerusalem,  Israel 

B.  Zaltzman,  Ben-Gurion  University  of  Negev,  Israel 
I.  Gotz,  Technische  Universitat  Munchen,  Germany 
A.  Novick-Cohen,  Technion  -  IIT,  Israel 

P.  Rybka,  Technische  Universitat  Munchen,  Germany 
W.  Merz,  Technische  Universitat  Munchen,  Germany 
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Arkady  Vilenkin 

Hebrew  University  of  Jerusalem,  Israel 

Effect  of  surface  dynamics  on  phase  transformation  and 

MASS  TRANSFER 

Three  different  problems,  combining  surface  dynamics  with  interface  motion  or  external  mass 
flux  through  the  interface  are  discussed.  The  first  describes  the  motion  of  the  phase  boundary 
in  two-dimensional  system  with  free  surfaces  corresponding  to  grooves  around  a  triple  junction. 
The  related  free  boundary  problem  consists  of  three  partial  derivative  equations  coupled  through 
boundary  conditions.  In  the  long  time  limit  the  problem  has  depending  on  input  data  either 
two  stable  traveling  wave  solutions,  or  a  life  time  of  the  system  is  limited.  The  second  problem 
concerns  the  motion  of  grain  boundary  in  a  two-dimensional  bicrystal  with  anisotropic  free  surfaces. 
Dufference  in  surface  energy  of  neighboring  grains  yields  grain  boundary  motion,  resulting  in 
decreasing  of  systems  energy.  Mathematical  structure  of  this  problem  is  similar  to  that  of  the 
previous  problem.  In  the  long  time  limit  is  studied.  In  this  limit,  depending  on  the  input  data, 
the  solution  is  either  a  stable  traveling  wave  or  does  not  exist  at  all.  A  peculiar  feature  of  this 
traveling  wave  solution  is  a  related  ”  negative  differential  resistance”  of  the  grain  boundary  velocity 
dependence  on  the  driving  force.  The  third  problem  describes  the  motion  of  the  free  surface  of  a  two 
dimensional  bicrystal  placed  between  two  rigid  walls.  The  respective  mathematical  formulation 
consists  of  a  fourth  order  partial  differential  equation  for  surface  dynamics  with  two  free  boundaries 
and  a  restriction  on  the  maximum  value  of  the  surface  shape  function  (the  requirement  of  non¬ 
overlap  of  the  two  crystals).  The  problem  possess  a  steady-state  solution  whenever  the  control 
parameter  (a  specified  mass  flux  along  the  grain  boundary  due  to  an  external  force)  is  below  some 
critical  value.  For  larger  fluxes  the  solution  oscillates  in  time,  with  each  oscillation  generating  a 
void,  which  propagates  along  the  grain  boundary  from  the  surface  into  the  bicrystal.  This  change 
of  morphology  decreases  the  effective  value  of  the  flux  imposed  on  the  free  surface. 


Boris  Zaltzman 

Ben-Gurion  University  of  the  Negev,  Israel 

Asymptotic  solutions  to  the  free  boundary  problem  arising 

IN  THE  grains  DYNAMICS 

The  boundary  value  problem,  arising  in  the  modelling  of  the  motion  of  grain  boundary  in  a 
two-dimensional  bicrystal,  has  been  studied.  The  studied  problem  is  a  free  boundary  problem  for 
a  fourth-order  partial  differential  equation.  The  well-posedness  of  this  free  boundary  problem  has 
been  studied  and  the  short  time  asymptotic  behavior  of  the  solution  analysed. 


I.  Gotz 

Technical  University  of  Munich,  Germany 

’’Undercooling”  in  type  i  superconductors 

We  consider  a  one-dimaensional  model  of  phase  changes  from  ’’normal”  to  superconducting 
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(see  [1],  [2]). 


ut  -  Uxx  =0,  0  <  a;  <  s{t), 
vt  -  Vxx  =0,  0  <  a:  <  s{t), 

Ux  =  —s'u  on  a;  =  s{t), 

Vx  —  —s'v  on  a;  =  s{t), 
u^+v^=H‘1  on  a;  =  s(t), 

Ux  =  Vx  =  on  X  —  0; 
u(x,0)  =  uo(x),  v(x,0)  —  vo(x), 

here  H  =  (u,  v)  is  the  magnetic  field  and  the  graph  of  the  function  s  separates  the  superconducting 
region  {a;  >  s{t)}  from  the  normal  one  {0  <  x  <  s(t)}.  A.  Friedman  and  B.  Hu  [2  ]  studied  this 
problem  under  the  assumption,  that  the  ’’normal”  region  is  continuously  increasing  in  time,  i.e.  the 
superconducting  region  shrinks.  We  are  interested  in  the  opposite  case,  when  the  superconducting 
region  is  increasing  in  time.  This  problem  is  closely  related  to  the  usual  one-phase  scalar  Stefan 
problem  with  undercooling  (consider  u  =  0).  That  problem  has  been  studied  in  the  paper  [3].  We 
have  proved,  that  under  some  critical  undercooling  a  blow-up  can  develope  and  the  free  boundary 
becomes  discontinuous.  The  question  was:  whether  the  similar  phenomena  are  usual  for  the  whole 
problem  (l)-(7).  Our  answer  is  :  NO.  Such  blow-ups  are  rather  an  exeption  for  our  vector  valued 
problem  (l)-(7).  More  exactly:  if  the  function  uq  is  more  than  some  positive  constant,  and  some 
smoothness  assumptions  are  satisfied,  then  the  function  s  is  Holder  continuous  and  piece-wise 
infinite  times  differentiable.  In  order  to  get  a  solution  we  take  a  kinetic  regularization: 

+v'^  =  +  es'  on  X  =  s{t), 

instead  of  (5).  Then  we  obtain  some  independent  on  e  estimates  for  the  functions  u,  v,  w  = 
u'^  +  v"^,  (j)  =  arctan{vlu)  and  tend  e  to  zero. 
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[2]  A.  Friedman  and  B.  Hu,  ”A  free  boundary  problem  arising  in  superconductor  modeling”, 
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Chaim  Charach  (Ben-Gurion  Univ.  of  the  Negev)  and  Paul  C.  Fife  (Univ.  of  Utah)  ’’Thermo¬ 
dynamically  consistent  phase-field  models  of  solidification  processes” 


Amy  Novick- Cohen 

Technion  -  IIT,  Israel 

Order-Disorder,  Phase  Separation,  and  Wetting 

A  system  of  Cahn-Hilliard/Allen-Cahn  equations  are  considered  which  model  simultaneous 
order /disorder  and  phase  separation.  It  is  demonstrated  that  this  system  exhibits  the  possibility 
of  triple  junction  motion  in  which  the  interface  between  the  two  ordered  phases  is  completely 
wetted  by  the  disordered  phase. 
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Piotr  Rybka 

Technische  Universitdt  Munchen,  Germany 

Asymptotics  for  an  equation  related  to  martensitic  phase 

TRANSITIONS 

We  study  asymptotic  behavior  of 

(1)  0  =  divcr(V'ii)  +  Aut  —  S^A^u, 

where  u  :  Q  C  IB?'  — >■  IR,  tr(^)  =  W’CC)  and  W  has  several  local  minima.  This  equation  is  the 
quasi-steady  approximation  of  the  viscoelasticity  with  capillarity 

(2)  putt  =  div(7(V'u)  -I-  Aut  —  5^A^u, 

Setting  p  =  0  corresponds  to  neglecting  the  inertial  effects  and  it  is  justifiable  for  slow  processes. 

Equation  (1)  arises  in  the  studies  of  van  der  Waals  fluids  and  phase  transitions  in  solids. 
Its  applicability  to  models  of  moving  phase  boundaries  has  been  justified  by  Slemrod  and  more 
recently  by  Abeyaratne  -  Knowles  in  the  context  of  phase  transitions  in  solids.  We  adopt  special 
boundary  conditions,  such  that  the  influence  of  the  boundary  on  the  overall  behavior  is  little.  We 
set  Cl  =  [0,  u)]  X  [0,  L]  and 


(3)  Ux2  =  Aux2  =0  at  X2  =  0,  L, 

We  define  the  energy  functional  to  be 

E{u)  =  f  [W(Vri(,r;)) 

Ja 

where  W{Fi,F2)  =  (f>{Fi)  +  ^F^,  and  </>  €  C°°{IR)  is  such  that,  (^(t)  >  0,  (p{t)  =  0  iff  \t\  =  1,  and 
Our  main  result  is 

Theorem.  Let  us  suppose  that  p  is  large,  W  satisfies  the  conditions  above  and  in  addition 
(j)  is  real  analytic  with  polynomial  growth.  If  uq  £  W^’^{Cl)  is  anti-periodic  satisfying  the  first 
of  boundary  conditions  (3),  then  there  exists  a  unique  strong  solution  of  (1)  which  satisfies  (3). 
Moreover,  if  in  addition  its  energy  E{uo)  is  sufSciently  small,  then  u{t)  converges  in  W^’^{Cl)  to 
a  minimizer  of  E. 

Let  us  point  out  that  because  of  the  boundary  condition  we  have  a  circle  of  minimizers.  The 
proof  of  this  Theorem  starts  with  a  study  of  critical  points  of  E.  We  show  that  global  minimizers 
do  not  depend  on  X2.  This  circle  of  minimizers  is  separated  from  other  critical  points  of  E  in 
the  topology  of  W‘^’^{Cl),  p  >  2.  It  turns  out  that  S'^E/Su^  has  a  two-dimensional  kernel  at  the 
minimizers.  Thus,  standard  tools  of  the  theory  of  dynamical  systems  are  not  applicable  here.  So, 
we  resort  to  methods  based  on  analyticity  of  W  and  relying  on  a  Lojasiewicz’  type  inequality. 
Joint  work  with  Karl-Heinz  Hoffmann. 


u  is  anti-periodic  in  xi. 

Y\^u{x)\'^]dx, 
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Wilhelm  Merz 

Technische  Universitdt  Munchen,  Germany 

Construction  of  an  asymptotic  model  for  the  oxidation 

PROCESS  OF  SILICON 

The  simulation  of  thermal  oxidation,  an  important  process  step  in  integrated  circuit  production, 
requires  the  computation  of  a  change  of  layer  structures  with  time.  At  a  process  temperature 
which  typically  lies  between  700-1200® C,  the  silicon  wafers  are  exposed  for  30-500  minutes  to 
a  gas  atmosphere  containing  oxygen  (O2).  Oxygen  enters  the  silicon-dioxide  layer  (Si02)  from 
outside,  diffuses  through  it  towards  the  silicon  (Si)  and  reacts  there  to  produce  new  silicon-dioxide 
according  to  Si  -|-  O2  — t  Si02. 

As  silicon  is  consumed,  the  boundary  layer  between  silicon  and  silicon-dioxide  moves  further  into 
the  silicon.  The  freshly  created  oxide  requires  more  than  twice  the  volume  of  the  consumed  silicon, 
for  which  reason  the  surface  of  the  oxide  layer  grows  outwards.  The  reaction  rate  and  the  volume 
expansion  are  usually  non-uniform  along  the  reaction  interface  and  this  results  in  non-homogeneous 
deformations  of  the  whole  layer  system. 

There  exists  a  variety  of  mathematical  models  describing  this  process.  We  consider  a  model,  where 
the  two  phases,  silicon  and  silicon-dioxide,  are  separated  by  a  moving  mixed  zone  of  thickness  0{e), 
£  being  a  small  parameter.  We  use  the  techniques  of  formal  asymptotic  analysis  to  derive  a  model 
with  a  sharp  reaction  front  between  the  two  phases.  The  resulting  free  boundary  problem  agrees 
with  other  existing  oxidation  models  on  the  basis  of  a  sharp  interface.  In  addition  we  give  a  proof 
for  the  existence,  uniqueness  and  qualitative  behaviour  of  the  solution  to  the  system  of  ODEs 
obtained  to  leading  order  of  the  interfacial  expansion. 
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Robert  Almgren 

University  of  Chicago,  USA 

Second-order  phase  field  asymptotics  for  unequal  conductivities 

The  phase  field  method  constructs  a  pair  of  stiff  reaction-diffusion  PDEs  whose  solutions  have 
a  thin  internal  layer.  The  PDE  system  is  constructed  so  that  as  a  parameter  in  the  equations 
tends  to  zero,  the  layer  approaches  a  sharp  interface,  satisfying  the  desired  boundary  conditions  at 
leading  order,  and  coupling  to  the  diffusion  equation  in  the  bulk.  In  numerical  computations,  it  is 
very  difficult  to  make  the  interface  thin  enough  to  achieve  quantitative  agreement  with  the  sharp- 
interface  model.  Recently,  Karma  and  Rappel  have  obtained  accurate  results  by  taking  account  of 
temperature  variations  within  the  interfacial  layer  and  by  slightly  modifying  the  kinetic  coefficient 
in  the  PDE.  We  show  how  their  construction  is  in  effect  a  second-order  asymptotic  analysis,  and 
we  extend  it  to  the  case  when  the  thermal  conductivities  are  different  in  solid  and  liquid,  as  for  the 
very  important  problem  of  alloy  solidification.  We  hope  that  second-order  asymptotic  analysis, 
and  identification  of  the  leading-order  error  terms,  will  soon  be  considered  to  be  as  fundamental 
as  are  their  analogs  for  numerical  discretization  schemes. 


Zhimirg  Chen 

Academia  Sinica,  China 

Numerical  methods  for  two  phase  continuous  casting  problems 

The  continuous  casting  problem  is  a  mathematical  model  describing  the  solidification  of  a 
material  being  cast  continuously  with  a  prescribed  velocity.  We  will  propose  a  finite  element 
scheme  which  is  based  on  an  implicit  discretization  of  the  diffusion  and  an  explicit  approximation  of 
the  convection.  For  the  problems  with  nonlinear  fiux  boundary  conditions  we  are  able  to  show  the 
converegence  of  the  method.  The  analysis  is  based  on  a  new  apriori  estimate  of  fractional  derivative 
in  time  of  the  discrete  temperature.  For  the  problems  with  linear  flux  boundary  conditions  we  are 
able  to  derive  an  error  estimate  of  the  same  convergence  order  as  that  analysis  is  based  on  a  new 
sharp  boundary  estimate  for  the  Green  operator  under  the  condition  that  the  free  boundary  does 
not  touch  on  the  inflow  boundary. 


G.  Georgiou 

University  of  Cyprus,  Cyprus 

Converged  solutions  of  the  planar  and  axisymmetric 

EXTRUDATE  SWELL  PROBLEMS 

We  solve  the  planar  and  axisymmetric  Newtonian  extrudate  swell  problems  using  standard  and 
singular  finite  elements.  We  study  the  convergence  of  the  two  methods  with  mesh  refinement,  for 
various  values  of  the  Reynolds  and  capillary  numbers.  The  singular  finite  elements  are  superior  to 
the  standard  finite  elements  only  when  coarse  or  moderately  refined  meshes  are  used,  the  Reynolds 
number  is  low  and  the  surface  tension  is  not  high.  The  standard  finite  elements  perform  better 
as  the  surface  tension  and/or  the  Reynolds  number  are  increased,  which  implies  that  the  effect  of 
the  stress  singularity  at  the  exit  of  the  die  becomes  less  important. 

Joint  work  with  Andreas  G.  Boudouvis  (  National  Technical  University  of  Athens). 
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Marianne  K.  Korten 

University  of  Buenos  Aires,  Argentina 

A  STRUCTURE  THEOREM  FOR  SOME  FREE  BOUNDARY  PROBLEMS  FOR 

THE  HEAT  EQUATION 

In  one  space  dimension  and  for  a  given  initial  datum  ui{x)  G  C^,  (say  such  that  Uf(x)  >  1  in 
some  interval)  the  equation  ut  =  A(u  -  1)+  can  be  thought  of  as  describing  the  energy  per  unit 
volume  in  a  Stefan-type  problem  where  the  latent  heat  of  the  phase  change  is  given  by  (l-u/(2;))  +  . 

Given  a  solution  in  the  sense  of  distributions  0  <  u  G  x  (0,T))  to  this  equation,  we 

prove  that  the  free  bounbary  is  n-rectifiable,  i.  e.,  but  a  null  set  the  union  of  images  of  imbedded 
manifolds  of  dimension  n  of  IR"  x  (0,T).  On  the  free  boundary  the  ’’classical”  Stefan  condition 
is  shown  to  hold  pointwise. 

We  show  the  same  result  for  the  combustion  model. 


A.  Myslinski 

System  Research  Institute,  Poland 

Fictitious  domain  approach  for  solving  shape  optimization 

PROBLEMS 

The  paper  deals  with  the  numerical  solution  of  a  shape  optimization  problem  of  an  elastic 
body  in  unilateral  contact  with  a  rigid  foundation.  The  equilibrium  state  of  this  contact  problem 
is  described  by  an  elliptic  variational  inequality  of  the  second  order. 

The  shape  optimization  problem  for  the  elastic  bodies  in  contact  consists  in  finding,  in  a 
contact  region,  such  shape  of  the  boundary  of  the  domain  occupied  by  the  body  that  the  normal 
contact  stress  is  minimized.  It  is  assumed  that  the  volume  of  the  domain  occupied  by  the  body 
is  constant.  Moreover  the  function  describing  the  boundary  of  the  domain  occupied  by  the  body 
and  its  derivative  are  bounded. 

Shape  optimization  of  elastic  or  hyperelastic  contact  problems  was  considered,  among  others 
in  [r5,r6,r7,r9,rll].  The  existence  of  optimal  solutions  was  investigated  in  [r9].  The  neccessary 
optimality  conditions  were  formulated  in  [rll]  using  the  material  derivative  approach  as  well  as  in 
[r5,r9]  using  the  optimal  control  theory  approach.  The  convergence  of  the  finite  element  approx¬ 
imation  was  investigated  in  [r6,r9].  Numerical  results  are  reported  in  [r6,r7,r9].  In  these  papers 
the  classical  approach  to  numerical  solving  of  optimal  shape  design  problems  based  on  boundary 
or  domain  variations  methods  [r4,r8,rl0]  was  used.  In  this  approach  the  state  problem  is  solved 
many  times  on  the  domain  changing  during  the  computation.  The  boundary  or  domain  variation 
methods  require  calculation  of  a  new  triangularization  of  the  optimized  domain,  updating  the 
stiffness  matrix  and  load  vector  at  each  iteration  of  the  numerical  algorithm.  Since  the  optimized 
domain  usually  has  complicated  geometrical  structure  the  whole  computational  process  is  tedious, 
time  consuming  and  expensive.  To  overcome  this  difficulty,  in  response  to  growing  number  of  in¬ 
dustrial  applications  of  the  optimal  shape  design  problems,  fixed  domain  methods  [r8]  for  solving 
these  problems  are  being  developed.  Fixed  domain  methods  are  based  on  using  fictitious  domain 
method  [r2,r3].  The  fictitious  domain  method  for  solving  the  state  systems  described  by  partial 
differential  equations  consists  in  transforming  the  original  state  system  defined  in  the  complicated 
geometry  domain  into  a  new  system  defined  in  a  given  fixed  simply  geometry  domain  containing 
the  original  domain  with  the  same  differential  operator  [r2,r3,r4].  This  method  allows  to  use  fairly 
structured  meshes  on  a  simple  geometry  domain  containig  the  actual  one  and  fast  solvers. 

Fixed  domain  methods  for  solving  shape  optimal  design  problems  were  considered  in  [r4,r8,rl0]. 
In  [r4]  the  original  shape  optimal  design  problem  for  the  Laplace  state  equation  with  Dirichlet 
boundary  condition  was  transformed  into  equivalent  one  using  a  new  control  variable  on  the  right- 
hand  side  of  the  state  equation  defined  in  the  fixed  domain  containing  the  original  one.  The 
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convergence  of  the  finite  element  approximation  is  shown  as  well  as  a  numerical  result  obtained  by 
using  nonsmooth  optimization  method  is  presented  in  [r4].  In  [r8]  the  original  shape  optimal  design 
problem  is  replaced  by  a  constrained  control  problem  with  either  boundary  control  or  distributed 
control.  The  state  system  is  defined  on  a  fixed  larger  domain  than  the  original  one.  The  existence 
of  optimal  controls  and  optimality  conditions  are  investigated  in  [r8].  The  numerical  algorithm 
based  on  the’  fixed  domain  method  for  shape  optimization  of  flow  problems  is  proposed  in  [rlO]. 

This  paper  deals  with  the  numerical  solving  of  a  shape  optimization  problem  for  linear  contact 
problem  using  the  fictitious  domain  approach.  To  author  knowledge  this  approach  for  contact 
problems  has  not  been  investigated  yet  and  the  present  paper  is  continuation  of  the  author  paper 
[r6] .  In  the  paper  we  shall  formulate  2D  contact  problem  in  variable  domain.  The  boundary  of  the 
domain  is  variable  subject  to  optimization.  The  cost  functional  approximating  the  normal  contact 
stress  is  introduced  [r6,r7].  The  fictitious  domain  formulation  to  this  optimization  problem  is 
given  using  the  boundary  control  technique  [r3,r9]  and  the  equivalence  to  the  classical  formulation 
is  shown.  The  piecewise  linear  finite  element  approximation  for  displacements  and  constant  for 
Lagrange  multipliers  is  proposed  and  its  convergence  is  shown.  The  design  sensitivity  analysis 
for  the  discrete  shape  optimization  problem  is  performed.  The  Lagrangian  multiplier  method  for 
solving  the  discretized  state  system  is  proposed  [r7].  The  descent  gradient  method  is  used  for 
solving  the  optimization  problem.  Numerical  examples  are  provided  and  discussed. 
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P.I.  Plotnikov 

Lavrentyev  Institite  of  Hydrodynamics,  Russia 

Singular  limits  of  solutions  to  Cahn-Hilliard  equation 

AND  AN  EVOLUTION  OF  MICROSTRUCTURES 
We  deal  with  the  singular  limits,  as  small  parameters  (r,  e)  0,  to  boundary-valued  problem 

(Ai9 - )t  =  'dxx,  rut  =  e^Uxx  -u^-'du  (0  <  a;  <  1,  0  <  t  <  T), 

2 

u(0,t)  =  =  '*9(0,  t)  =  i5(l,t)  =  0, 

u{x,0)  =  uo{x),  d{x,0)  =  doix),  (0  <  X  <  1). 

which  may  be  regarded  as  linearization  of  Penrose-Fife  model  near  the  critical  temperature.  It  is 
supposed  that  r  =  e“,  a  G  (1,2)  and  initial  data  satisfy  the  following  conditions: 

e|lwoi||L2(o,i)  +  11^*0, i9o||c(o,i)  <  c,  uo,i9o  C  0^(0, 12). 

Under  these  assumptions  the  problem  has  the  unique  smooth  solutions  (i9e,U£).  The  basis  idea  is 
that  a  description  of  the  passage  to  limit  can  be  down  in  terms  of  slow  variables:  the  temperature 
de  and  the  adiabatic  invariant  le  =  ‘i.xiul  -t-  i^e)^  -  The  main  result  is  the  following: 

Theorem  Let  {ue,'ds)  — >■  {u,'d)  weakly  in  L2{Q)-  Then: 

a) {u£,'de)  (u,t9)  strongly  in  L^iQ)-  and  the  inclusion 

€{I,d:d<0,0<i<d^}LS{I  =  ^2}. 

holds  almost  everywhere  in  Q  =  (0, 1)  x  (0,T) 

b)  Whenever  continuous  even  function  g{u), 

/OQ 

K{z,  /,  •d)g{z)dz 

-OO 

where  K  is  some  algebraic  function. 


B.  Stoth 

Universitdt  Bonn,  Germany 

A  DEGENERATE  PARABOLIC  ELLIPTIC  SYSTEM  AND  A  FREE 
BOUNDARY  PROBLEM 

In  a  Type-II  superconductor  the  magnetic  field  penetrates  the  superconducting  body  through 
the  formation  of  vortices.  Because  the  number  of  vortices  is  large  it  seems  feasible  to  model  their 
evolution  by  an  averaged  problem,  known  as  the  mean  field  model  of  superconductivity.  In  this 
talk  I  will  discuss  a  two  dimensional  reduction,  assuming  all  vortices  to  be  perpendicular  to  a  given 
direction.  Since  both  the  magnetic  field  H  and  the  averaged  vorticity  w  are  curl  free,  they  may 
be  represented  via  a  scalar  magnetic  potential  q  and  a  scalar  stream  function  ijx,  respectively.  I 
will  present  existence  and  uniqueness  results  of  solutions  (i/’)9)  of  fh®  resulting  degenerate  elliptic 
parabolic  system 

^t-|VV'l(aV-(||$y)+9-V')  =0  inflx]0,r[, 

-Ag -f- xnq  =  Xnf/'  iniR2x]0,T[, 
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by  means  of  viscosity  and  weak  solutions.  In  addition  I  will  relate  (■0,q)  to  solutions  of 

the  mean  field  equations.  Finally  I  will  discuss  the  regularity  properties  of  the  corresponding 
stationary  free  boundary  problem  and  present  special  solutions  to  it. 

This  is  joint  work  with  C.  M.  Elliott  and  R.  Schatzle. 


D.  Tarzia 

Univ.  Austral,  Argentina 

Numerical  analysis  of  the  steady-state  two-phase 
Stefan-Signorini  problem 

We  consider  a  steady-state  heat  conduction  problem  in  a  multidimensional  bounded  domain 
with  regular  boundary.  We  assume,  without  loss  of  generality,  that  the  melting  temperature  is 
zero  degree  centigrade.  We  consider  a  source  term  g  in  the  domain,  a  convective  condition  on 
a  part  B1  of  the  boundary,  a  positive  heat  flux  q  on  a  part  B2  and  its  remain  part  B3  of  the 
boundary  we  have  a  Signorini  type  condition. 

We  study  the  numerical  analysis  and  the  corresponding  error  bounds  between  the  continuous 
and  the  discrete  condition  in  order  to  obtain  a  two-phase  Stefan-Signorini  problem  (that  is,  the 
solution  is  a  function  of  non-constant  sign  in  the  whole  domain)  through  the  numerical  analysis 
of  the  variational  inequalities  by  using  the  finite  element  method. 


Alfred  Wagner 

Institut  fur  Mathematik  -  Koeln,  Germany 

The  instationary  Bernoulli  problem 

We  consider  an  inviscid,  incompressible  and  irrotational  fluid  in  a  bounded  two  dimensional 
domain  O  x  [0,T].  Imbedded  in  the  fluid  (described  by  the  potential  u  )  we  imagine  a  gas  bubble 
of  prescribed  volume  and  pressure.  The  interface  gas  -  fluid  is  assumed  to  be  of  constant  pressure 
i.e.  we  do  not  include  surface  tension.  The  shape  of  the  bubble  is  described  as  a  level  set  of  a 
function  to  be  determined. 

We  prove  the  existence  of  a  global  generalised  solution  in  classical  spaces  using  the  technique  of 
generalised  sub  -  and  supersolutions.  The  solution  may  not  be  unique  due  to  possible  fattening 
effects  of  the  interface. 
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A.  Yu.  Zubarev 

Ural  State  University,  Russia 

Corrugation  of  ferronematic  free  surface 

Ferronematics  are  colloidal  suspensions  of  ferromagnetic  particles  in  a  liquid  crystal.  Unusual 
and  outlooking  for  high  technologies  properties  of  these  compounds  attract  considerable  attention 
of  researchers.  In  well-known  papers  the  behaviour  of  ferronematics  in  the  volumes  with  the 
fixed  boundaries  has  been  studied.  The  aim  of  this  study  is  an  analysis  of  the  behaviour  of 
ferronematic  free  boundary  in  magnetic  field  and  in  its  absence.  Two  types  of  ferronematics  have 
been  considered,  the  first  one  being  compensated  (this  means  that  magnetization  in  the  absence 
of  the  field  equals  zero)  and  the  second  one  being  uncompensated  (with  the  magnetization  being 
constant  in  its  absolute  value).  The  situations  have  been  studied  when  ferronematic  director  is 
normal  and  parallel  to  its  boundary.  The  critical  values  of  normal  and  tangential  field  above  which 
the  plate  surface  becomes  unstable  has  been  derived.  In  the  supercritical  fields  the  hexagonal  or 
quadratic  lattice  of  valleys  and  ridges  is  developed;  the  growth  of  deformations  is  limited  by 
capillary  forces  and  Fkank  tensions. 

If  the  field  applied  is  normal  to  ferronematic  director  and  exceeds  some  critical  value,  then  a 
deformation  of  a  free  surface  reveals  in  conjunction  with  the  space  division  of  ferroparticles  the 
moments  of  which  are  directed  opposite  to  each  other.  In  the  most  real  cases  the  corrugation  of  a 
free  surface  is  a  soft  phase  transition;however,  in  certain  cases  it  may  occur  as  a  rigid  transition. 
It  should  be  noted  that  the  deformations  of  paramagnetic  or  dielectric  liquids  free  surface  in  static 
tangential  field  has  not  been  yet  discussed  in  literature. 

If  magnetization  of  uncompensated  (ferromagnetic)  ferronematics  is  normal  to  its  surface  and 
exceeds  critical  value,  then  a  free  surface  of  such  system  becomes  corrugative  and  hexagonal 
lattice  of  valleys  and  ridges  occurs  on  this  surface  even  in  the  absence  of  the  field.  The  critical 
values  of  external  field  and  ferronematic  magnetization  have  been  determined  from  the  balance  of 
capillary,  magnetic  and  gravitational  forces  as  well  as  Frank  tensions  in  a  liquid  crystal  matrix. 
The  amplitude  of  the  free  surface  corrugation  is  calculated  as  well  as  its  period  and  the  type  of 
the  occuring  lattice.  Depending  on  the  system  parameters,  the  corrugation  of  ferromagnetic  liquid 
crystal  free  surface  may  appear  as  a  soft  or  rigid  phase  transition. 
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Toyohiko  Aiki 

Gifu  University,  Japan 

Large-time  behavior  of  two-phase  Stefan  problems  for 
NONLINEAR  PARABOLIC  EQUATIONS 


We  consider  the  following  two-phase  Stefan  problems:  Find  a  function  u  =  u{t,x)  on  Q(T)  = 
(0,  T)  X  (0, 1),  0  <  T  <  oo,  and  a  curve  x  —  i{t),  0  <  ^  <  1  on  [0,  T]  satisfying  that 


p{u)t  -  a{ux)x  -f  ^  +  9{u)  = 


fo  in 

fl  in  Q^p{T), 

f_{t,x)  £  (3{u{t,x))  for  a.e.  {t,x)  £  Q{T), 

Qf\T)  =  {{t,x);0  <t<T,0<x<  1(1)}, 
g  W  (T)  =  {(t,  x);  0  <  f  <  T,  i{t)  <x<l}, 
u{t,  l{t))  =0  for  0  <  t  <  T, 
it)  = -a{ux){t,t{t)-)  +  iov  a..e.  t  £  [0,T], 

a{ux){t,0+)  £  dbQ{u{t,0))  for  a.e.  t  £  [0,T], 

-  a{ux){t,l~)  £  db{{u{t,l))  for  a.e.  t  G  [0,T], 
u{0,x)  =  uo{x)  for  X  £  [0,  l]i 

^(0)  =  4, 


where  p  :  R  -i-  R  and  a  ■.  R—^  R  are  continuous  increasing  functions;  /3  is  a  maximal  monotone 
graph  on  R-,  g  :  R  R  is  a,  Lipschitz  continuous  function;  fi{i  =  0, 1)  is  a  given  function  on 
(0,  oo)  X  (0, 1);  =  0, 1)  is  a  proper  l.s.c.  convex  function  on  R  for  each  t  >  0  and  dbt  denotes 

its  subdifferential  in  i?;  uo  is  a  given  initial  function  and  £o  is  a  number  with  0  <  ^0  <  1- 

In  this  paper  we  discuss  global  existence  in  time  of  solution,  global  attractor  and  w-limit  set 
(cf.  [1]).  Recently,  Ito,  Kenmochi  and  Yamazaki  [2]  obtained  the  interesting  results  concerned 
with  the  attractor  for  non-autonoumous  problems  for  the  evolution  equations  in  a  Hilbert  space. 
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Asymptotic  behavior  of  solutions  to  phase  field  equations 

WITH  CONSTRAINTS  UNDER  NONLINEAR  DYNAMIC  BOUNDARY 

CONDITIONS 

We  consider  a  nonlinear  system  of  the  following  form:  Find  functions  u  :  [0,  -foo)  — >■  L‘^{Cl),w  : 
[0,  -f-oo)  and  v  :  [0,  -boo)  L^(r)  satisfying  that 

ut  +wt  —  Au  =  f{t,  x)  in  Q  (0,  -t-oo)  x  fl, 

owt  —  K  Aw  +  (i{w)  -b  g{w)  B  u  in  Q, 
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(1) 

(2) 


(3) 


u  —  V  a.e.  on  S  :=  (0,  +oo)  X  F, 


— -f  c— + /i(a)  =  0  on  S, 
an  at 

(4) 

dw 

—  =  0  on  S, 
an 

(5) 

u{0,x)  =  uo{x),  w{0,x)  —  wo{x)  for  x  £  Q, 

(6) 

v{0,x)  =  vo{x)  for  a;  G  F, 

(7) 

where  is  a  bounded  domain  in  with  smooth  boundary  F;  u,  k  and  c  are  positive  constants; 

/3  is  a  maximal  monotone  graph  in  R  x  R',  g  :  R  R  and  h  :  R  R  are  Lipshitz  continuous 
functions;  /  is  a  given  function  on  Q;  uq,  wq  and  vo  are  initial  functions.  ^  denotes  the  outward 
normal  derivative  on  F.  We  denote  by  CP=CP(uo,  raoi 'i^o)  the  above  system  (1)  (7). 

The  boundary  condition  (4)  for  the  unknown  function  u  is  called  the  nonlinear  dynamic  bound¬ 
ary  condition,  which  includes  the  time  derivative  of  the  unknown  function  on  the  fixed  boundary 
and  is  regard  as  a  mathematical  model  for  a  feedback  control  device  on  the  fixed  boundary. 

In  the  present  talk  we  discuss  the  existence  and  uniqueness  of  the  solution  to  CP,  noting  that 
the  condition,  c  >  0,  is  required  in  order  to  proVe  the  uniqueness  by  using  the  usual  method  for 
the  evolution  equation  on  the  dual  space  of 


G.Bayada(^),  M.Chambat^^^  and  C.Vazquez^^) 

^INSA,  France  -  “^Univ.  Lyon  I,  France  -  ^Univ.  of  La  Coruna,  Spain 
Numerical  analysis  of  a  cavitation  free  boundary  problem 

WITH  FLUX  IMPOSED  BOUNDARY  CONDITIONS 

In  this  work,  an  upwinding  semidiscretized  scheme  based  on  characteristics  method  is  ana¬ 
lyzed  when  it  is  applied  to  a  free  boundary  problem.  The  free  boundary  problem  is  issued  from 
a  cavitation  model  in  lubrication  and  it  involves  especific  fiux  imposed  boundary  conditions.  The 
boundary  conditions  correspond  to  an  axially  lubricated  journal  bearing  device  and  the  math¬ 
ematical  analysis  is  the  subject  of  the  previous  work.  In  this  paper  we  prove  that  the  specific 
semidiscretization  method  leads  to  a  well  posed  family  of  semidiscretized  problems  that  converges 
to  the  departure  continuous  problem. 


Piotr  Biler 

University  of  Wroclaw,  Poland 

Global  in  time  versus  blowing  up  solutions  to 

CHEMOTAXIS  SYSTEMS  WITH  FUNCTIONAL 
COEFFICIENTS 

We  study  global  in  time  solvability  for  initial  data  of  arbitrary  size  to  some  simplified  versions 
of  the  Keller-Segel  model  of  chemotaxis 


lit  =  V  •  (Vu  —  uV(l){v)), 
evt  =  Au  —  V  +  u. 
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This  parabolic  {e  =  1)  or  parabolic-elliptic  (e  =  0)  system  is  considered  in  a  bounded  smooth 
domain  C  under  the  homogeneous  Neumann  conditions 


du 


dv 

dv 


=  0  on  X  (0,  oo), 


and  suitable  initial  conditions.  The  system  (*)  is  a  particular  example  of  systems  describing  the 
movement  of  microorganisms  sensitive  to  the  gradient  of  a  chemical  secreted  by  themselves.  Here 
u  =  u{x,t)  >  0  and  v  =  v{x,t)  >  0  represent  the  densities  of  the  population  and  of  the  chemical, 
X  £  Q,  t  >  0.  The  function  (p  is  called  the  sensitivity  function,  so  that  V(j){v)  is  the  velocity  of  the 
chemotactic  movement. 

Recently,  much  attention  has  been  attracted  by  the  basic  model  with  (f>{v)  =  v  and  e  =  0.  In 
particular,  the  questions  of  global  in  time  continuation  of  solutions  with  small  initial  data,  versus 
finite  time  blow-up  of  large  ones  have  been  considered.  Radially  symmetric  solutions  that  explode 
after  a  concentration  of  mass  have  been  constructed  by  M.  A.  Herrero  and  J.  J.  L.  Velazquez.  Nu¬ 
merical  evidence  of  blow-up  at  the  boundary  has  been  reported  by  K.  Gajewski  and  K.  Zacharias. 
Quantitative  conditions  (in  terms  of  concentration  of  initial  data  measured  by  e.g.  their  Morrey 
space  norms)  have  been  given  by  the  author  for  a  related  model  of  gravitational  interaction  of 
particles.  An  extensive  list  of  problems  under  current  study  include  characterization  of  blowing 
up  solutions,  a  description  of  free  boundary  of  the  region  of  spreading  of  u  and  a  comparison  of 
chemotaxis  systems  with  different  sensitivity  functions  (j>. 

Here  we  generalize  some  results  of  T.  Nagai  and  T.  Senba  who  considered  radially  symmetric 
solutions  of  the  system  (*),£  =  0,  in  balls,  give  extensions  to  the  nonradial  case,  and  an  extension 
to  a  parabolic  system  of  chemotaxis  with  e  =  1. 

T.  Nagai  and  T.  Senba  studied  the  sensitivity  functions  4i{v)  =  with  the 

constants  x  >  0  and  0  <  p  <  1.  Their  results  included  (for  the  both  functions  (p)  the  global  in 
time  existence  of  radial  solutions  in  the  disk  in  the  plane  {N  =  2),  and  for  (p{v)  —  with 

X  <  2/(A^  ~  2)  —  in  the  balls  in  ,  N  >  3.  Moreover,  for  N  >  3  and  either  logarithmic  (p 
with  X  >  2N/{N  —  2)  or  (p{v)  —  xr'^i  there  exist  radial  solutions  that  blow  up  in  a  finite  time.  In 
certain  sense,  (p{v)  =  xlogu  is  a  critical  nonlinearity  for  the  existence  of  global  solutions  to  (*). 

THEOREM  1.  Suppose  that  N  =  2  and  the  sensitivity  function  (p  is  strictly  sublinear;  0  < 
(p' {v)  — >  0  when  v  oo,  as  well  as  v<p'(v)  is  increasing.  Then  each  radially  symmetric  solution 
to  (*)  with  e  =  0  and  a  bounded  density  u{x,0)  can  be  continued  to  the  global  in  time  one. 


THEOREM  2.  If  N  >  2,  (p{v)  =  xlog'^  <^'^-<1  x  <  2/N’,  then  the  system  (*)  with  e  =  0  and 
arbitrary  initial  data  u(x,0)  £  L^(fl)  (for  some  q  >  N/2)  has  global  in  time  solutions. 

THEOREM  3.  If  N  —  2  and  (p{v)  —  xlogu  with  x  S  !>  then  for  all  initial  data  the  system  (*) 
with  e  =  1,  u{x,0)  £  and  Vi;(a;,0)  £  has  global  in  time  solutions. 


A.C.  Briozzo,  M.F.  Natale  and  D.A.  Tarzia 

Univ.  Austral,  Argentina 

Determination  of  unknown  thermal  coefficients  through 

PARABOLIC  FREE  BOUNDARY  PROBLEMS 

We  determine  unknown  thermal  coefficients  of  a  semi-infinite  material  through  a  phase-change 
process  with  an  overspecified  condition  on  the  fixed  faceby  considering  two  parabolic  free  boundary 
problems.  In  both  problems  we  obtain  formulae  for  the  unknown  thermal  coefficients  and,  the 
necessary  and  sufficient  conditions  for  the  existence  of  an  explicit  solution  of  the  corresponding 
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free-moving  boundary  problem.  In  the  first  part  we  consider  a  material  of  Storm’s  type  with  non¬ 
constant  thermal  coefficient.  Here,  we  study  a  free  boundary  and  a  moving  boundary  problem 
for  the  heat  equation  (Stefan-like  problem)  for  9  cases.  In  the  second  part,  we  consider  a  free 
boundary  problem  for  a  non-linear  heat  conduction  equation  with  a  convective  term  for  14  cases. 


Nikolaus  Bubner 

Berlin,  Germany 

Modelling  and  control  of  martensitic  phase  transitions 
IN  shape  memory  alloys 

Shape  memory  alloys  (SMA;  e.g.  CuZnAl,  NiTi,  ...  )  have  received  increasing  attention  in 
recent  years  partly  due  to  the  wide  range  of  their  technical  applic  ability  in  actuators  and  adaptive 
structures.  They  show  a  noticeable  change  in  their  mec  hanical  behaviour,  i.e.  stress-strain  relation 
in  different  temperature  ranges:  elastic  at  high  tempera  tures  and  pseudo-  or  quasiplastic  at  low 
temperatures.  An  alloy  can  be  permanently  de  formed  to  up  to  10%  without  fracture  and  still 
it  recovers  its  old  shape  under  heating  or  co  oling.  This  is  due  to  first  order  structural  phase 
transitions  between  different  equilibrium  configurations  of  the  metallic  1  attice,  named  austenite 
and  martensite.  Austenite  is  the  undeformed  crystal  lattice  which  is  stable  at  high  temperatures. 
By  deforming  the  lattice,  one  obtains  24  crystallographically  equivalent  variants.  These  martensite 
variants  prevail  at  low  temperatures. 

In  order  to  exploit  these  phenomena  for  industrial  applications  a  satisfactory  m  odel  is  needed. 
We  develop  a  thermodynamically  consistent  one-dimensional  Landau-Ginzburg  model  which  is 
able  to  describe  all  relevant  effects  [1,2].  For  instance,  we  investigate  a  so-called  deformation-driven 
experiment  (DDE):  a  thin  rod  of  a  SMA  is  fixed  on  one  side  and  pushed  and  pulled  on  the  other 
side  in  the  course  of  time  by  an  elongation  m.  To  demonstrate  the  validity  of  our  model  we  will 
compare  the  numerical  results  to  experiments  performed  by  I.  Muller  and  his  co  -workers  [3,4]. 
Summarizing,  for  a  DDE  we  have  the  following  system  {Cl  :=  {0,1),  Q  :=  Cl  X  (0; T)): 

(1)  putt  -  {^{d  -  0i)ux  -  Pul +aul)x  +  Suxxxx  =  0,  in  Q> 

(2)  Ce  Ot  —  tvBxx  y  0  Ux  Uxt  ~  9 {^^l) :  in  G , 

u{0,t)  —  Uxx{0,i)  —■  Uxx{lil)  “  d?  u{l,t)  —  TTl{t), 
ex{o,t)  =  o,  -K9x{i,t)  =  K{e{i,t)-er{t)),  vtG[o,T], 
u(a:,  0)  =  uo(a:),  'at(a;,  0)  =  ui  (i),  d{x,0)  =  9o{x),  Va;  G  12. 

(1)  and  (2)  represent  the  balance  laws  of  momentum  and  energy,  respectively.  It  denotes:  p 
-  mass  density,  u  -  displacement,  9  -  absolute  te  mperature,  Ux  -  strain,  Ce  -  specific  heat,  k  - 
heat  conductivity,  g  -  density  of  heat  sources  or  sinks,  I  -  length  of  the  rod,  R,  -  heat  exchange 
coefficient,  6r  -  outside  temperat  ure.  The  couple  stress  leads  to  the  Ginzburg-term  6  ■  Uxxxx, 
and  a,  P,  7,  and  d  are  material  constants  to  be  determined  for  each  specimen. 

We  discuss  existence  and  uniqueness  results  for  (1)  (see  [5]),  we  outline  the  determination  of 
the  physical  parameters,  and  we  present  our  nume  rical  simulations  [2].  The  main  features  are: 

1)  Displacement,  strain  and  temperature  evolution  show  the  experimentally  observed  phenomena: 
phase  transitions  between  variants  o  f  martensite  or  between  austenite  and  one  variant,  respec¬ 
tively,  changes  of  temperatures  and  latent  heats,  nucleation  processes,  propagation  of  p  hase 
boundaries. 

2)  Size  and  shape  of  our  simulated  hysteresis  loops  are  independent  of  dicret  ization.  The  coeffi¬ 
cient  of  the  term  modelling  the  interfacial  energy,  i.e.  6,  det  ermines  the  size  of  the  loop  and 
leads  to  quantitative  agreement  with  experiments.  This  corresponds  to  the  model  of  [3]. 

3)  The  model  is  even  capable  of  producing  interior  hysteresis  loops  which  are  also  found  in  the 
experiments. 
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In  order  to  control  the  system,  we  will  consider  a  weak  formulation  of  (1)  {WF).  The  aim  is 
to  achieve,  possibly  isothermally,  a  prescribed  distribution  of  the  phases.  Therefore,  it  is  natural 
to  consider  a  cost  functional  involving  the  order  parameter  e  —  Ux  and  6  (or  rather,  the  stress 
a  =  7(0  —  9i)ux  —  +  OiU%),  as  well  as  the  natural  control  variables  m,  6^,  and  g.  For  instance, 

one  could  investigate  the  case  when 

J{m,g,er)  =  ai||  a  -  ct  +  “2II  ^  ^  |li2(Q) 

+  «3llhi||^2(o,T)  +  «4||  5  lli2(Q)  +  as II  fi'r  lli2(o,T)’ 

where  ai,  i  =  1, . . . ,  5,  are  non-negative  constants,  and  where  a  and  6  denote  the  desired  stress 
and  temperature  distribution  duri  ng  the  evolution  of  the  process,  respectively.  Then  our  control 
problem  reads 

(CP)  Minimize  J(m,g,  0r),  subject  to  WF  and  {m,g,6r)  G  Uad- 

Here,  Uad  denotes  the  set  of  admissible  controls,  and  is  some  nonempty,  convex,  bounded,  and 
closed  set. 

We  prove  the  Frechet  differentiability  of  the  solution  operator  and  we  derive  the  necessary 
conditions  of  optimality  [5].  Of  course,  also  other  cost  functionals  are  conceiv  able  in  actual 
applications. 

Furthermore,  investigating  a  more  realistic  case,  we  include  the  following  state  constraints  for 
stress  and  temperature  [6],  i.e. 

{6,cr)  e  C,  with  C  given  by 

C  :=  I  {9,a)  e  C{Q)  X  C{Q)^ci  <  d{x,t)  <  C2,  C3  <  a{x,t)  <  c^,  V  (x,  t)  6  Q  |. 

Again,  the  necessary  optimality  conditions  are  given. 

As  an  application  we  are  going  to  dampen  out  the  vibrations  of  a  thin  rod  by  hea  ting  and 
cooling.  We  want  to  emphasize  the  fact  that  all  these  kind  of  problems  can  be  regarded  as 
subproblems  for  actuators  and  adaptive  structures. 
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POINTWISE  AND  VISCOSITY  SOLUTIONS  FOR  THE  LIMIT  OF  A  TWO  PHASE 
PARABOLIC  SINGULAR  PERTURBATION  PROBLEM 

In  [CLWl],  [CLW2]  we  deal  with  the  problem  of  studying  the  uniform  properties  and  the  limit, 
as  e  — >■  0,  of  solutions  u^{x,t)  of  the  equation 

(Ps)  Au^ -ui  = 

where  £  >  0,  de  >  0,  Be{s)  =  support  d=[0,l]  and  J  f3{s)ds  =  M. 

We  consider  a  family  of  solutions  to  in  a  domain  X>  C  1  which  are  uniformly  bounded 
in  L°°  norm  in  T>.  We  get  uniform  estimates  that  allow  us  to  pass  to  the  limit  as  £  ^  0,  we  find 
properties  of  the  limit  function  u  in  general  situations,  and  we  prove  that  u  is  a  solution  to  the 
free  boundary  problem 


Au  -ut  =  0  in  \  d{u  >  0} 
u  =  0  ,  {u+f-{u-f^2M  on  Vr\d{u>  0}. 

Here  u  is  the  inward  unit  spatial  normal  to  the  free  boundary  VDdiu  >  0},  it+  =  max(u,  0)  and 
u~  =  max(— u,  0). 

On  one  hand,  we  prove  that  every  limit  function  u  is  a  solution  to  P  in  a  very  weak  sense. 
Moreover,  at  any  free  boundary  point  (xo,to)  there  holds  that  if  limsup^j,^^).^^^,^  IVu  |  <  7, 

then  limsup(^_t)_^(3.p_tQ)  \Vu+\  <  \/2M  +  72. 

On  the  other  hand,  we  prove  that,  in  the 

two  phase  case,  the  free  boundary  condition  is  satisfied  in  a  pointwise  sense  at  any  “regular” 
free  boundary  point  (xo,fo)'  We  point  out  that  we  only  make  assumptions  on  u  at  (xo,to)-  More 
precisely,  if  {xo,to)  €  P  n  d{u  >  0}  is  such  that  there  exists  an  inward  unit  spatial  normal  to  the 
free  boundary  V  fl  d{u  >  0}  in  a  parabolic  measure  theoretic  sense  and  such  that  {u  =  0}  has 
zero  “parabolic  density”  at  {xo,to)  then,  there  exist  a  >  0,7  >  0  such  that 

u{x,  t)  —  a{x  -  XQ,v)'^  -  7(x  -  xo,v)~  +  o(|a;  —  xo|  +  |t  —  to|  ^  ) 

with  —  7^  =  2M.  Here  v  is  the  inward  unit  spatial  normal  to  the  free  boundary  at  (xo,to). 

Finally  we  prove  that,  when  {u  =  0}°  =  0,  the  limit  function  u  is  a  viscosity  solution  to  the 
free  boundary  problem  P.  By  a  viscosity  solution,  we  mean  a  weak  solution  of  the  free  boundary 
problem  in  the  sense  that  it  is  a  continuous  function  for  which  comparison  principles  with  classical 
supersolutions  and  subsolutions  hold. 
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Localization  properties  for  a  quasilinear  degenerated 
SYSTEM  IN  Semiconductors  Theory 

In  this  joint  work  with  J.I.  Diaz  and  A.  Jiingel,  we  consider  the  following  problem:  let  u 
(electrons  density),  v  (holes  density)  and  w  (electric  potential)  be  any  solution  of  the  system 

{ut  —  div  {V(p{u)  —  b{u)V'w)  =  F{u,  v) ' 
vt  —  div  (V(p(v)  +  b(v)Vw)  =  F(u,  v)  ►  . 

—Aw  =  V  —  u  +  C. 

Find  conditions  on  (p  and  F  that  give  rise  to  the  existence  of  free  boundaries.  We  use  a  local 
energy  method  combined  with  the  consideration  of  the  characteristics  defined  by  Vta  to  get  the 
following  results:If  (p(s)  :=  s^,  with  m  >  1  then  if  F(s,<t)  satisfies 

[^(3,  cr)(s’^  +cr"^)  <  c  ,]  then  u,v  have  the  property  oi  finite  speed  of  propaga¬ 
tion.  Moreover,  if  F(s,cr)  satisfies  [—F(s,  (7)(s”‘ +0'’^)  >  c  ,  with  0< 

p  <  ^  —  1]  and  the  global  energy  is  bounded  then,  even  if  the  initial  data  are  strictly  positive,  a 
dead  core  arises. 
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Domain  decomposition  methods  for  the  stefan  problem  in  non 

HOMOGENEUS  MEDIA:  SOME  INDUSTRIAL  APPLICATIONS 

Nonoverlapping  domain  decomposition  methods  are  very  adecuated  to  solve  problems  where 
certain  ”  difficult”  phenomena  occurs  only  in  a  specific  subdomain,  whereas  in  the  rest  of  the  domain 
a  standar  equation  applies.  This  can  be  the  case  of  the  Stefan  problem  in  non  homogenous  media. 
A  lineal  model  problem  is  formulated  as  an  optimization  problem  with  constrains,  and  so  reduced  to 
the  localization  of  a  saddle  point  of  a  Lagrangian  or  Augmented  Lagrangian  funtionals.  This  idea  is 
due  to  Glowinski  and  Le  Tallec  ( 1990)  but  is  fully  developed  in  this  work.  We  present  a  wide  family 
of  nonoverlapping  decomposition  methods  and  study  their  convergence  properties  for  the  linear 
case.  We  also  present  some  numerical  results  for  the  Stefan  problem  in  non  homogenous  media 
related  to  the  following  industrial  applications:  a)The  aluminium  solidification  in  a  cast.  b)The 
deodorization  of  vegetable  oils  which  involves  condensation  and  freezing  of  a  multicomponent  vapor 
mixture  inside  shell  and  tube  heat  exchanger. 
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Growth  of  the  visual  image  in  electrophotography 

Friedman  and  Velazquez  [1]  have  modelled  the  development  of  the  visual 
image  in  electrophotography  as  a  time-dependent,  free  boundary  problem.  They  proved  results 
on  existence  and  uniqueness  of  a  solution  to  a  two-dimensional  version  of  the  problem.  We  present 
a  numerical  solution  to  the  same  version  of  the  problem  and  find  the  shape  of  the  toner  layer  (i.  e. 
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the  visual  image)  as  a  function  of  time.  At  the  initial  instant,  an  exact  expression  can  be  derived 
for  the  electric  field  and,  hence,  for  the  growth  rate  of  the  toner  layer.  We  use  that  expression 
to  start  the  computation.  We  also  derive  an  integral  equation  formulation  of  the  free  boundary 
problem,  and  use  that  formulation  to  find  the  shape  of  the  toner  layer  at  later  times. 

[1]  Friedman,  A.  and  Velazquez,  J.L.:  A  time  dependent  free  boundary  problem  modeling  the 
visual  image  in  electrophotography,  Arch.  Rational  Mech.  Anal.  123  (1993),  259-303. 
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Two  LIMITING  REGIMES  OF  THE  DIRECTIONAL  SOLIDIFICATION 
OF  BINARY  MELTS  WITH  THE  MUSHY  REGION 

The  concentrational  supercooling  of  a  binary  melt  before  crystallization  front  results  in  the 
origination  of  mushy  region.  An  analysis  of  the  structure  of  the  mushy  region,  which  includes 
the  liquid,  solid  particles,  and  dendrites  extending  from  the  bulk  solid  surface,  is  suggested.  The 
analytical  solution  of  the  quasi-equilibrium  two-phase  zone  problem  is  obtained  on  the  basis  of  the 
small  parameter  method;  the  temperature,  concentration  of  solute  inside  the  region,  and  solid  phase 
volume  fraction  are  found.  If  the  zone  thickness  is  much  smaller  than  the  pertinent  linear  scales 
of  the  process  under  consideration,  the  zone  itself  may  be  substituted  by  a  discontinuity  interface. 
A  new  formulation  of  the  corresponding  mathematical  problem  with  an  unknown  boundary  is 
suggested,  which  modifies  the  known  Stefan  problem  by  means  of  replacing  the  mass  balance 
condition  at  a  smooth  interface  with  a  new  one  stemming  from  the  analysis  of  the  zone  structure 
and  enables  to  derive  physically  transparent  conclusions. 

When  heat  supply  to  a  binary  melt  at  its  outer  boundary  is  absent,  macrokinetics  of  solidifica¬ 
tion  are  governed  by  progressive  cooling  of  the  melt  rather  than  by  withdrawal  of  an  admixture. 
Because  the  former  process  is  much  faster  than  the  latter  one,  a  metastable  region  of  concentra¬ 
tional  supercooling  arises  and  has  time  to  expand  before  the  ingot  surface  noticeably  changes  its 
position.  This  leads  to  bulk  crystallization  within  such  a  region  prior  to  solidification  at  the  surface. 
If  crystalline  grains  are  forming  and  growing  rapidly  enough,  the  two-phase  region  is  practically 
at  the  state  of  thermodynamical  equilibrium.  A  simple  self-similar  solution  of  releavant  equations 
is  obtained  which  approximately  decsribes  unidirectional  process. 
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Stability  for  solutions  of  Allen-Cahn  equation 

We  consider  the  asymptotic  behaviour  of  Allen-Cahn  equation  with  constraint  in  one¬ 
dimensional  space: 

tnt  —  (l?nx|^~^n;x)x  +  C  +  p(u^)  =  0  in  Q, 

^  ^  in  Q, 

subject  to  initial  and  boundary  conditions; 


Wx{t,  ±L)  =0  for  t  >  0, 
w{0,x)  =  woix)  for  any  x  G  {—L,L). 
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Here  L  >  0  is  a  positive  number;  p  is  a  constant  with  2  <  p  <  +oo;  g{w)  is  a  smooth  function 
of  w]  a*]  subdifferential  of  the  indicator  function  g..]  of  the  compact  interval 

[a*,  a*]  C  R. 

For  this  model,  we  consider  the  following  second-order  elliptic  problem  as  steady-state 
problem; 

-(!«;?> r2ro~)x+^°°+p(u;~)=0  in  (-L,L), 

in  (-L,L), 
wf{±L)  =  0. 

In  this  talk,  we  discuss: 

(i)  structure  of  the  set  of  stationary  solutions  u;°°. 

(ii)  convergence  of  w{t)  to  a  stationary  solution  w°°  as  t  — >■  -foo. 

(hi)  stability  of  stationary  solutions  u;°°. 
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Viscosity  solutions  and  regularity  of  the  free  boundary 

FOR  THE  LIMIT  OF  AN  ELLIPTIC  TWO  PHASE  SINGULAR 
PERTURBATION  PROBLEM 

In  this  work  we  are  concerned  with  the  following  problem,  of  interest  in  combustion;  Study  the 
limit  as  e  ->  0,  of  solutions  (x)  to  the  equation 

(Be) 

where  e  >  0,  Peis)  =  >  0  in  (0, 1),  yS  =  0  otherwise,  and  f  P(s)ds  =  M. 

Given  a  family  of  uniformly  bounded  solutions  to  Be  in  a  domain  C  one  expects  that 
the  limit  function  u  be  a  solution  to  the  free  boundary  problem: 

Art  =  0  in  f2  \  d{u  >  0}, 

^  ^  u  =  0,  =  2M  on  Q  n  >  0}, 

where  u  is  the  inward  unit  normal  to  the  free  boundary  H  n  >  0}.  In  fact,  in  this  work  we 
prove  the  following  result  s  for  the  limit  function  u: 

1)  Let  xo  G  d{u  >  0}  be  such  that  there  exists  an  inward  unit  normal  o  to  d{u  >  0}  at  xq  in 
the  measure  theoretic  sense.  If  either  the  set  {u  <  0}  has  positive  density  or  u"*"  is  nondegenerate 
at  xo  then  u  satisfies  the  free  boundary  condition  at  xq.  I.e.,  there  exist  both  ut  and  uj,  and 

—  [uZ)"^  =  2M'  at  XQ. 

2)  Under  suitable  assumptions,  u  is  a  solution  to  J5  in  a  viscosity  sense.  By  a  viscosity  solution 
we  mean  a  continuous  function,  harmonic  away  from  the  free  boundary,  which  satisfies  the  free 
boundary  condition  in  terms  of  proper  asymptotic  developments.  In  particular,  u  is  a  viscosity 
solution  if  {u  =  0}°  =  0  or  if  is 

in  terms  of  proper  asymptotic  developments.  In  particular,  u  is  a  viscosity  solution  if  {u  = 
0}°  =  0  or  if  is  nondegenerate  on  d{u  >  0}. 

3)  We  prove,  in  addition,  the  following  regularity  results: 

3.i)  Under  suitable  assumptions  on  the  limit  function  u,  there  is  a  subset  of  the  free  boundary 
d{u  >  0}  which  is  locally  a  surface.  Thus,  on  this  portion  of  d{u  >  0},  the  free  boundary 
condition  is  satisfied  in  the  classical  sense.  Moreover,  this  smooth  subset  is  open  and  dense  in 
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d{u  >  0}  and  the  remainder  of  the  free  boundary  has  {N  —  l)-dimensional  Hausdorff  measure 
zero. 

3.ii)  If  u~  is  locally  uniformly  nondegenerate  on  d{u  >  0},  then  the  free  boundary  is  locally  a 
surface.  Hence,  there  are  no  singularities  and  u  is  a  classical  solution  to  E.  We  point  out  that 
the  loca  1  uniform  nondegeneracy  of  u~  on  d{u  >  0}  is  a  necessary  condition  for  the  regularity 
of  the  free  boundary  in  the  strictly  two  phase  case,  even  if  we  only  require  the  free  boundary  to 
be  locally  Lipschitz  continuous. 

We  finally  want  to  remark  that  there  are  limit  functions  u  which  do  not  satisfy  the  free  boundary 
condition  in  the  classical  sense  on  any  portion  of  d{u  >  0}  (for  instance,  u  —  axf  +  ax^  with 
0  <  a  <  y/2M). 


C.  Lederman^,  J.  L.  Vazquez^  and  N.  Wolanski^ 

^  Univ.  de  Buenos  Aires,  Argentina  -  ^  Univ.  Autonoma  de  Madrid,  Spain 

Uniqueness  of  solution  to  a  free  boundary  problem  from  combustion 

We  report  the  results  in  [LVW]  where  we  investigate  the  uniqueness  and  agreement  between 
different  kinds  of  solutions  for  a  free  boundary  problem  in  heat  propagation  that  in  classical 
terms  is  formulated  as  follows:  to  find  a  continuous  function  u{x,t)  >  0,  defined  in  a  domain 
V  C  X  (0,  T)  and  such  that 


ut  —  Au  in  V  r\  {u  >  0}. 

Besides  we  assume  that  the  interior  boundary  of  the  positivity  set,  V  fi  d{u  >  0},  so-called  free 
boundary,  is  a  regular  hypersurface  on  which  the  following  conditions  are  satisfied 

a  =  0,  —du/dv  =  C 

{v  the  outward  unit  spatial  normal).  In  addition,  initial  data  are  specified  ,  as  well  as  either 
Dirichlet  or  Neumann  data  on  the  external  boundary  of  V  H  {u  >  0}.  This  problem  arises  in 
combustion  theory  as  a  limit  situation  in  the  propagation  of  premixed  flames  (high  activation 
energy  limit). 

This  free  boundary  problem  admits  classical  solutions  only  for  good  data  in  the  small.  Several 
generalized  concepts  of  solution  are  proposed,  among  them  the  concepts  of  limit  solution  and 
viscosity  solution.  We  investigate  conditions  under  which  the  three  concepts  agree  and  produce  a 
unique  solution. 

A  limit  solution  is  a  function  which  is  obtained  as  the  limit  (as  e  — >■  0)  of  solutions  to 


Au^  -  uf  = 

where  0  <  /3e  <  supp  Be  =  [0,  e]  and  f  /3e(s)  ds  =  ^. 

A  viscosity  solution  is  a  continuous  function  which  satisfies  comparison  principles  with  classical 
supersolutions  and  subsolutions  to  the  free  boundary  problem. 

Our  results  can  be  summarized  as  saying  that,  under  suitable  assumptions  on  the  domain, 
the  reaction  function  Be  and  on  the  initial  and  boundary  data,  if  a  classical  solution  to  the  free 
boundary  problem  exists  in  a  certain  time  int  erval,  then  it  is  at  the  same  time  the  unique  classical 
solution,  the  unique  limit  solution  and  also  the  unique  viscosity  solution  in  that  time  interval. 
For  related  works  see  [CV],  [CLW],  [GHV]. 
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A.  Mamourian 

Univ.  of  Tehran,  Iran 

On  a  non-uniformly  elliptic  system  boundary 

VALUE  PROBLEM 

Linear  boundary  values  problems  for  a  degenerate  nonlinear  system  of  the  Lavrentiev  type 
equations,  in  the  case  when  the  index  corresponding  to  the  boundary  values  problems  are  negative, 
has  been  investigated.  Making  use  of  functiontheoretic  methods  in  the  theory  of  partial  differential 
equations,  existence,  uniqueness  and  stability  of  the  solution  of  boundary  problems  in  complex 
form,  for  non-uniformly  ellipticity  of  the  equation  in  the  Sobolev  space  are  introduced. 
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Karol  Mikula 

Slovak  Technical  University,  Slovakia 

Solution  of  nonline arly  curvature  driven 

EVOLUTION  OF  PLANE  CURVES 

We  investigate  the  evolution  of  closed  smooth  planar  curves  T  :  R/Z  B?  satisfying  the 
geometrical  equation 

(1)  v  =  0{k) 

where  v  is  the  normal  velocity  of  evolving  curve  and  /3  is  a  smooth  function  of  the  curvature  k.  As 
a  typical  example  one  can  consider  a  function  P{k)  =  k^  where  m  >  0.  Geometrical  equations  of 
type  (1)  are  used  for  describing  various  phenomena  in  physics,  material  sciences,  computer  vision, 
robotics  and  artificial  intelligence.  There  are  two  main  fields  in  which  the  evolution  of  plane  curves 
play  an  important  role: 

a)  the  multiscale  analysis  of  images  and  shapes  closely  related  to  signal  smoothing,  edge  detection, 
segmentation  and  image  representation; 

b)  the  Stefan  problem  with  surface  tension  and  related  interface  motion  models. 

We  suggest  a  computational  method  for  solving  (1).  The  aim  is  to  represent  the  equation  by  a 
so-called  intrinsic  heat  equation 


(2) 


dx  _  d^x 
dt  ds^ 


where  s*  is  an  appropriate  parametrization  of  a  curve  corresponding  to  the  function  /3.  Then, 
we  make  use  of  the  ’Eulerian  transformation’  of  the  intrinsic  heat  equation  into  a  degenerate 
evolutionary  partial  differential  equation  with  spatial  variable  being  independent  on  time  and 
varying  on  a  fixed  interval.  This  equation  is  a  generalization  of  the  corresponding  one  studied  by 
Dziuk  for  classical  curve  shortening  flow  problem.  We  prove  some  a-priori  estimates  of  a  smooth 
solution  which,  in  particular,  imply  the  curve  shortening  property  of  the  governing  equation. 
The  same  property  is  proven  for  time  semidiscretization.  The  proposed  numerical  scheme  is 
carefully  tested  by  various  examples  of  the  nonlinear  curvature  driven  evolution  (1).  Numerical 
results  are  compared  with  the  exact  selfsimilar  solutions  and  also  with  previous  results  obtained 
by  a  conceptually  different  method  introduced  by  Mikula  and  Kacur  (which  can  be  used  only  in 
convex  case) .  We  apply  the  method  for  computing  the  evolution  of  nonconvex,  selfintersecting  and 
nonrectifiable  curves  as  well. 

This  is  a  joint  work  with  Daniel  Sevcovic. 


Andreas  Miinch 

Technical  University  of  Munich,  Germany 

The  DESTABILISING  EFFECT  OF  SURFACE  NON-UNIFORMITIES  ON 
THIN  FLUID  FILMS  FLOWING  DOWN  AN  INCLINATION 

In  recent  studies  it  could  be  shown  that  linear  stability  analysis  for  a  thin  film  flowing  down  an 
inclined  plane  predicts  stabilisation  below  a  certain  critical  inclination  angle  a*  depending  on  the 
static  contact  angle  6s-  Experiments,  however,  have  shown  Anger  formation  for  inclination  angles 
significantly  below  a*  in  situations  where  9s  is  small. 

Since  every  realistic  plane  surface  exhibits  non-uniformities  at  a  microscopic  scale,  the  possi¬ 
bility  that  these  may  serve  as  an  additional  source  of  instability  is  considered  in  this  work. 
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We  model  the  non-uniformities  as  a  perturbation  of  the  slip  length,  a  parameter  which  enters 
in  a  singular  fashion  into  the  governing  lubrication  equation  for  the  thin  film. 

On  the  basis  of  our  model  we  present  numerical  evidence  that  these  perturbations  are  amplified 
to  macroscopic  changes  in  the  fluid  profile  and  the  contact  line,  supporting  the  proposition  that 
this  effect  should  play  an  important  role  explaining  the  mechanism  of  destabilisation  for  the  regime 
of  small  inclination  and  contact  angles. 


C.M.  Murea 

Univ.  of  Bucharest,  Romania 

Optimal  control  approach  for  the  fluid-structure  interaction 

PROBLEMS 

A  fluid-structure  interaction  problem  is  studied.  We  are  interested  by  the  displacement  of  the 
structure  and  by  the  velocity  and  the  pressure  of  the  fluid. 

The  contact  surface  between  fluid  and  structure  is  unknown  a  priori,  therefore  it  is  a  free 
boundary  like  problem. 

In  the  classical  approaches,  the  fluid’s  and  structure’s  equations  are  coupled  via  two  boundary 
conditions:  equality  of  the  fluid’s  and  structure’s  velocities  at  the  contact  surface  (which  is  a 
Dirichlet  like  boundary  condition)  and  equality  of  the  fluid’s  and  structure’s  forces  at  the  contact 
surface  (which  is  a  Neumann  like  boundary  condition). 

In  our  approach,  the  Dirichlet  like  contact  boundary  condition  will  be  relaxed  and  treated  by 
the  Least  Squares  Method. 

We  start  with  a  guess  for  the  contact  forces.  The  displacement  of  the  structure  can  be  computed. 
We  suppose  that  the  fluid’s  domain  is  completely  determined  by  the  displacement  of  the  structure. 
Knowing  the  actual  domain  of  the  fluid  ant  the  contact  forces,  we  can  compute  the  velocity  and 
the  pressure  of  the  fluid. 

In  this  way,  the  Neumann  like  contact  boundary  condition  is  trivially  accomplished. 

The  problem  is  to  find  the  contact  forces  such  that  the  Dirichlet  like  contact  boundary  condition 
holds.  It’s  a  exact  controllability  problem  with  Neumann  like  boundary  control  with  Dirichlet  like 
boundary  observation.  The  control  appears  also  in  the  coefficients  of  the  fluid’s  equations. 

In  order  to  obtain  some  existence  results,  this  exact  controllability  problem  will  be  transformed 
in  an  optimal  control  problem  using  the  Least  Squares  Method. 

This  mathematical  model  permits  to  solve  numerically  the  coupled  fluid-structure  problem  via 
partitioned  procedures  (i.e.  in  a  decoupled  way,  more  precisely  the  fluid’s  and  the  structure’s 
equations  are  solved  separately). 

The  aim  of  this  paper  is  to  present  the  existence  of  an  optimal  control  for  a  fluid-beam  inter¬ 
action  problem. 
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Fernando  Quiros  and  Juan  Luis  Vazquez 

Universidad  Autonoma  de  Madrid,  Spain 

Asymptotic  convergence  of  the  Stefan  problem  to 

Hele-Shaw 


A  Stefan  problem  is  a  mathematical  model  for  describing  the  melting  of  a  body  of  ice  in  contact 
with  a  region  of  water.  In  one-phase  Stefan  problems,  the  temperature  of  the  ice  is  supposed  to 
be  maintained  at  0°C.  Hence,  the  unknowns  are  the  temperature  distribution  in  the  water  and 
the  shape  of  the  interface  (free  boundary)  between  ice  and  water. 

The  Hele-Shaw  problem  is  a  model  describing  the  movement  of  a  viscous  fluid  confined  in  a 
narrow  cell  between  two  parallel  plates.  It  has  applications  in  the  plastics  industry  (injection 
moulding)  and  in  electromachining.  It  can  also  be  considered  as  the  zero  specific  heat  limit  of  the 
Stefan  problem,  being  in  that  sense  a  simplified  version  of  that  famous  problem. 

The  aim  of  this  work  (detailed  proofs  are  given  in  [2])  is  to  study  the  large-time  behaviour  of 
the  solution  to  the  initial-and-boundary-value  problem  in  an  exterior  domain,  Cl,  for  both  models, 
and  to  show  their  relationship.  We  prescribe  Dirichlet  boundary  data  g{x)  on  the  compact  set  5Q, 
which  are  non-negative,  non-trivial  {g  ^  0)  and  constant  in  time.  For  the  Stefan  problem  we  also 
specify  initial  data  6{x,  0)  =  do{x),  which  are  assumed  to  be  nonnegative  and  compactly  supported 
in  II.  For  the  Hele-Shaw  problem  we  need  only  specify  the  initial  support,  which  is  bounded.  The 
classical  formulation  of  these  problems  in  spacial  dimension  n  =  1  is  the  following: 

1-D  Stefan  problem.  To  find  9{x,t)  and  r{t)  such  that 


^  ae  _  d‘^6 

dt  dx"^ 
9{x,Qi)  =  9o{x) 

9{r{t),t)  =  0 
ox 

6>(0,t)  =  A 
^  r(0)  =  b. 


0  <  X  <  r{t),  t  >  0, 
0  <  X  <b, 

t  >  0, 
t  >  0, 
t  >  0, 


1-D  Hele-Shaw  problem.  To  find  p{x,t)  and  r{t)  such  that 


d'^p 

dx^ 


=  0 


p{r{t),t)  =  0 


^(r(t),f)  =  -ir'(t) 
ax 

p{0,t)  =  A 


1  r(0)  =  b. 


0  <  X  <  r{t),  t  >  0, 
t  >  0, 
t  >  0, 
t  >  0, 


We  will  show  that  there  is  a  marked  difference  in  the  asymptotic  behaviour  of  these  problems 
in  one  and  several  dimensions.  Thus,  though  in  one  space  dimension  there  is  convergence  to  a 
stationary  state  in  the  x/t^/^-scale  for  both  problems,  the  asymptotic  profile  is  not  the  same. 
The  situation  is  different  for  n  >  1.  Indeed,  there  is  no  self-similar  solution  which  explains 
simultaneously  the  large-time  behaviour  of  the  solutions  and  interfaces  of  our  problems.  The  actual 
behaviour  can  be  described  using  the  ideas  of  matched  asymptotic  expansions.  The  positivity  set 
expands  to  cover  in  finite  time  any  compact  subset  of  the  exterior  domain  Cl  and,  in  the  usual 
(x,  t)  coordinates,  the  solutions  of  both  problems  stabilize  to  the  solution  of  the  exterior  Dirichlet 
problem  for  the  Laplacian.  This  is  what  we  call  the  near  held  or  regular  limit.  On  the  other  hand, 
we  prove  that,  after  a  scaling  of  the  form 

y  =  V  =  if  n  >  3, 

y  =  x/TZ{t),  V  =  ulogTZ{t)  if  n  =  2, 
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where  'R.{t)  ~  C'(t/logt)^/^  as  t  oo,  and  u  represents  the  temperature  if  we  are  studying  the 
Stefan  problem  and  the  pressure  if  we  are  considering  the  Hele-Shaw  flow,  we  have  convergence 
in  both  problems  to  the  same  free  boundary  problem,  which  corresponds  to  a  singular,  radial, 
self-similar  solution  of  the  Hele-Shaw  problem.  This  is  what  we  call  the  far  Held  or  singular  limit. 
This  means  that  in  the  far  region,  close  to  the  free  boundary,  Hele-Shaw  represents  the  reduced 
problem  for  Stefan  (in  the  sense  of  Prandtl’s  theory  of  boundary  layers).  In  practice,  the  heat 
equation  of  the  Stefan  problem  Joses  the  inertial  term  as  t  ^  oo  because  this  term  becomes 
lower  order  with  respect  to  the  main  terms  at  the  asymptotic  level.  It  is  thus  a  case  of  asymptotic 
simpliHcation,  a  rather  usual  phenomenon  in  nonlinear  parabolic  equations.  The  analysis  gives  in 
particular  the  position  of  the  free  boundary,  whose  behaviour  is  given  by 

|a;|  ~  if  n  >  3, 

|a;|  ~  C{t/  \ogt)^^^  if  n  =  2. 

The  matching  between  both  regions  amounts  to  adjusting  the  constant  C  in  the  far  region. 

A  previous  study  of  the  asymptotic  behaviour  for  the  Stefan  problem  for  n  >  3  was  done  by 
Matano  in  [1].  He  proves  that  any  weak  solution  becomes  eventually  classical,  that  is,  the  interface 
between  the  ice  and  the  water  regions  is  sufficiently  smooth  for  all  large  t.  He  also  proves  that  the 
shape  of  the  free  boundary  (that  is,  the  interface)  tends  to  a  sphere  of  radius  (Ct)^^^  as  t  — >■  oo 
and  calculates  the  constant.  In  addition  to  that,  we  establish  the  asymptotic  developments  for  all 
n  >  2.  We  also  prove  that  the  shape  of  the  free  boundary  of  the  Hele-Shaw  problem  approaches 
a  sphere,  with  precisely  the  same  radius. 

In  conclusion,  there  is  a  strong  difference  between  the  one-dimensional  and  the  multidimensional 
cases:  for  n  =  1  solutions  of  the  Stefan  problem  do  not  eventually  become  solutions  of  Hele-Shaw, 
while  for  n  >  2  they  do.  There  is  a  explicit  selfimilar  solution  giving  the  asymptotic  behaviour 
of  Stefan  when  n  =  1  which  is  not  a  solution  of  Hele-Shaw.  For  n  >  2  we  need  two  matched 
developments  to  explain  the  asymptotic  behaviour  of  the  solutions  in  both  problems.  For  the 
record,  the  development  in  the  region  close  to  the  free  boundary  is  always  self-similar.  The 
influence  of  the  compactly  supported  initial  data  disappears  (to  first  order)  in  all  cases. 
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L.  Preziosi 

Politecnico  di  Torino,  Italy 

Non-Isothermal  injection  moulding  with  resin  cure, 

PREFORM  DEFORMABILITY,  AND  COUPLING 
BETWEEN  INTERFACED  DOMAINS 

The  paper  deduces  a  three-dimensional  model  based  on  the  theory  of  deformable  porous  media 
aimed  at  simulating  some  injection  moulding  processes  used  to  fabricate  composite  materials.  It 
works  under  non  isothermal  conditions,  includes  resin  cure  and  allows  the  solid  constituent  in 
both  the  dry  and  the  wet  region  to  deform  during  infiltration.  The  model  is  simplified  in  the  one¬ 
dimensional  case  by  performing  analytically  the  integrations  of  the  mechanical  equations  in  the 
uninfiltrated  region.  The  remaining  system  of  partial  differential  equations  in  the  two  interfaced 
and  time-dependent  domains  is  then  posed  with  the  proper  interface  and  boundary  conditions 
both  in  the  case  of  given  inflow  velocity  and  in  the  case  of  given  pressure  cycle.  After  writing  for 
numerical  reasons  the  problem  in  a  Lagrangian  formulation  fixed  on  the  solid  constituent  domain 
decomposition  techniques  are  used  for  the  simulation. 


A.  Shidfar  and  Karim  Ivaz 

Iran  University  of  Science  and  Technology,  Iran 


In  this  paper  we  consider  an  inverse  problem  of  heat  equation  with  free  boundary  condition. 
We  identify  necessary  initial  and  boundary  temperature  to  create  a  certain  free  boundary. 


George  G.  Tsypkin 

Moving  boundary  problems  of  geothermal  reservoir 

MODELING 

Coupled  processes  of  heat  and  multi-phase  mass  transfer  with  phase  transitions  in  porous  media 
are  of  practical  interest  in  geophysical  and  engineering  sciences.  A  large  class  of  these  problems 
are  Stefan  -  like  problems.  A  new  mathematical  models  of  water-ice-vapor  phase  transitions 
in  frozen  soils  [1],  gas  hydrate  dissociation  in  strata  [2],  and  phase  transition  water-steam  in 
geothermal  reservoirs  [3]  have  been  proposed.  These  models  are  generalizations  of  the  Stefan 
problem  and  take  into  account  appearance  of  two  unknown  moving  phase  transitions  boundaries 
and  (or)  extended  phase  transition  zones  (’’mushy  regions”).  The  models  include  equations  of 
mass  and  energy  conservation  laws,  the  generalized  Darcy  law  and  thermodynamic  relations.  The 
formulation  of  each  problem  contains  moving  interface  at  which  the  saturation  function  suffers 
a  discontinuity.  An  essential  part  of  the  mathematical  description  is  the  correct  specification  of 
boundary  conditions.  The  conservation  laws  may  be  written  in  term  of  ’’Jump”  at  phase  transition 
fronts.  We  assume  that  all  phase  at  the  fronts  are  locally  at  a  thermodynamic  equilibrium.  These 
relations  constitute  the  complete  system  of  boundary  conditions  at  the  moving  interface.  On 
the  whole  the  problems  are  nonlinear.  It  is  vital  to  emphasize  that  in  all  zones  it  is  possible  to 
apply  linear  approach  and  derive  the  self-similar  analytical  solutions  for  half  space  domain.  The 
solutions  in  all  regions  can  be  expressed  in  terms  of  linear  combinations  of  probability  integrals. 
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By  substituting  these  solutions  into  the  boundary  conditions  at  the  moving  fronts  we  obtain  the 
systems  of  transcendental  equations  for  unknown  functions,  which  are  solved  numerically.  For 
example,  the  problems  of  steam  (steam-water  mixture)  production  from  geothermal  reservoirs  are 
considered.  It  is  shown  that  the  introduction  of  an  evaporation  surface  separating  the  steam  and 
water  leads  to  superheating  of  the  water  in  a  zone  ahead  of  the  front.  This  contradiction  is  removed 
by  introducing  of  an  extended  phase  transition  zone  (’’mushy  region”)  between  the  single-phase 
zones.  Water  and  steam  in  this  region  are  locally  at  a  thermodynamic  equilibrium.  The  boiling 
front  separating  the  water  zone  and  the  two-phase  one  is  degenerated.  In  this  case  condition  of 
continuity  of  water  saturation  function  is  satisfied.  This  function  at  the  mixture-steam  front  varies 
discontinuously  from  the  unknown  value  to  zero. 
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Existence  of  periodic  solution  for  phase  transitions 

We  consider  a  solid-liquid  phase  transition  model. 

Let  9  be  the  absolute  temperature  and  w  be  the  order  parameter.  Then  our  model  is  a  couple 
of  system  of  nonlinear  second  order  parabolic  PDEs  (kinetic  equations  for  internal  energy  and  the 
order  parameter)  as  follows: 

[9  +  A(iu)]t  -  A  +  p9  =  f{t,x)  in  Q  (0,  -l-oo)  x  O, 

wt  -  kAw  -f  /3{w)  +  g{w)  +  ^  3  0  in  <3, 

U 

with  boundary  conditions 

■  |;(-i)+n„(-l)=M(.»).  1^=0  onS;:=(0,+oo)xr, 
and  To-periodic  conditions 


9{0,-)  =  9{To,-),  w{0,-)  =  w{To,-)  on  fl. 

Here  f2  is  a  bounded  domain  in  i?^(l  <  AT  <  3)  with  smooth  boundary  P  :=  5f2;  f3{-)  is  a  mciximal 
monotone  graph  in  Rx  R-,  A(-)  is  a  smooth  convex  function  on  R-,  g{-)  is  a  smooth  function  on  R] 
no,  K  and  p  are  positive  constants;  To  is  a  positive  constant  and  /,  h  are  given  data  which  are 
To-periodic  in  time. 

In  this  talk,  we  show  the  existence  of  periodic  solution  of  our  model. 
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The  simulation  of  instabilities  for  two- phase  flow 

BETWEEN  PARALLEL  PLATES 

The  hydrodynamical  instabilities  are  interesting  phenomena,  which  are  important  for  different 
technological  processes.  We  investigate  the  instabilitie  s  of  two-phase  flow  between  two  parallel 
plates.  The  mathematical  model  is  based  on  well  known  approximation  of  the  Hele-Shaw  flow. 
We  solve  the  problem  with  finite  difference  method.  For  describing  of  the  moving  free  boundary 
it  is  convenient  to  use  the  level-set  approach  [C]. 

The  level-set  method  needs  values  of  the  speed  at  grid-points  nearby  the  interface.  Because  the 
viscosity  and  the  pressure  have  jumps  on  the  interface,  it  is  an  open  problem,  how  to  approximate 
the  coupling  conditions  and  calculate  the  speed.  We  compare  two  different  approximations.  One  is 
based  on  the  smoothed  delta-function,  the  other  approximation  is  based  on  ideas  of  the  immersed 
interface  method  [B],[A].  The  proposed  methods  allow  the  simulation  of  instabilities  that  lead  to 
changes  of  the  topology  of  phases.  We  investigate  both  methods  and  present  numerical  results. 
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